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by 
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Submitted to the Department of Ocean Engineering and the Department of Mechanical 
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Master of Science in Naval Architecture and Marine Engineering 

and 

Master of Science in Mechanical Engineering 


ABSTRACT 

The thesis compares the analytical solution, two marine classification society design 
rules, and two design guides against experimental results for predicting the failure modes 
(general instability, axisymmetric buckling, and asymmetric collapse of the shell [lobar 
buckling]) and failure pressures of ring-stiffened cylinders 

The analytical solution is first summarized based on several sources. The design rules for 
the classification societies and the design guidance from two sources are then presented with 
brief explanations for each one. The design rules used are: American Bureau of Shipping ( Rules 
for Building and Classing Underwater Vehicles, Systems, and Hyperbaric Facilities, 1990) and 
Germanischer Lloyd ( Rules for Underwater Technology, 1988). The design guides used were 
Society of Naval Architects and Marine Engineers {Submersible Vehicle Systems Design, 1990) 
and Massachusetts Institute of Technology Course 13A Professional Summer Notes (MIT 13A 
Submarine Design Trends, 2001). 

The United States Navy Naval Sea Systems Command, Submarine Structural Integrity 
Division supplied experimental data for four cylinders that covered the failure modes and 
allowed comparison between experiment and design rules / guidance. 

The comparison of experimental to predicted data found that the design codes and design 
guides performed adequately in predicting axisymmetric yield and asymmetric buckling. The 
performance of the design codes and guides in predicting failure by general instability was 
unsatisfactory. For the experimental failures by general instability, the design codes and guides 
predicted significantly higher failure pressures than those experimentally determined; resulting in 
the design codes and guides actually predicting failure by axisymmetric yield in stead of general 
instability. These inconsistencies in the predictions of failure mode and pressures for general 
instability should be further explored to determine causes and corrections. 
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Chapter 1: Introduction 


The widespread use of stiffened cylinders in the marine industry has generated many 
studies into the stability and failure of these cylinders and methodology for failure prevention. 

Of primary concern to entities involved with the use of manned submersible vehicles is the 
design of ring-stiffened cylinders; this type of stiffened cylinder is used for significant portions 
of the pressure hull. Over a hundred years of theoretical and experimental research has led to a 
general understanding of the mechanics of failure for these cylinders. Based on this research, 
marine Classification Societies, such as the American Bureau of Shipping (ABS), the American 
Petroleum Institute (API), NORSOK and Germanischer Lloyd (GL) have promulgated design 
rules to provide guidelines on the design and building of stiffened cylinders for safe operation. 
Other design and analysis theories and guidance are available in texts such as those published by 
the Society of Naval Architects and Naval Engineers (SNAME). 

1.1 Failure Types 

There are three primary types of failure of ring-stiffened cylinders. They are 
axisymmetric yielding (AY) of the shell between stiffeners, asymmetric buckling of the shell 
between stiffeners (Lobar), and general instability of the shell and stiffeners (GI). Axisymmetric 
yield generally occurs when the shell is relatively heavy and the frames are closely spaced. 

Lobar buckling can occur when the shell is relatively light and the frames are strong and widely 
spaced. General Instability can occur when the cylinder is relatively long, the shell is thin and 
the frames are small. General Instability is very dependent upon eccentricities in the shell, which 
tend to lower the cylinder’s resistance to the General Instability mode. [1] 
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As analyzed here General Instability is presumed to occur in the elastic region of the 
stress - strain curve of the material. Cylinders also fail by inelastic General Instability, which 
occurs at significantly lower pressures than that of elastic General Instability. Failure by this 
mode is not addressed by the design rules. Other modes of failure also exist such as multi-wave 
instability, which is a sub-type of General Instability. It can occur in both the elastic and 
inelastic regions. Again the design rules do not address this failure mode. Several of the 
classification societies address stiffener tripping, which is the rotation of a stiffener away from 
perpendicular with the shell, however stiffener tripping is usually a precursor to general 
instability and is not a separate major failure mode. 

1.2 Concept Exploration 

For this thesis, an emphasis was placed on exploring how the various design rules 
predicted failure of cylinders that replicate modern submarine design (i.e. the shell was relatively 
thick compared to the diameter of the cylinder). This was facilitated by experimental failure data 
from the U.S. Navy’s Naval Sea Systems Command (NAVSEA) Submarine Structural Integrity 
Division. The analysis was limited to ring-stiffened cylinders to remain consistent with the 
primary concern of the thesis. The test cylinders that were chosen had failed in all three possible 
modes, allowing for comparison of the design rules in all modes of failure. . 

1.3 Analysis Techniques 

For the analyses, the design rules for the various classification societies and design 
guidance were programmed into MATHCAD™ for consistency of approach, ease of symbolic 
representation, and quickness of computations. The analytical solution was also programmed 
into MATHCAD™ for comparison. The experimental data from the four test cylinders were 
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then input into each computer code. The codes gave a failure pressure for each type of failure. 
The lowest calculated pressure was then considered the failure pressure and the corresponding 
mode was designated the failure mode. The failure modes and pressures were then compared 
against the experimental results with emphasis placed on agreement between codes and 
experiment on failure mode (first priority) and then failure pressure (second priority). If a 
different failure mode was predicted than that experimentally found for a particular cylinder, an 
analysis was performed for agreement between predicted failure pressures and also the closeness 
of the failure mode pressures. 

1.4 Design Rules and Guides Examined 

There were two classification society design rules examined: The American Bureau of 
Shipping (rules from Rules for Building and Classing Underwater Vehicles, Systems and 
Hyperbaric Facilities, 1990 Edition) (reference 2) and Germanischer Lloyd {Rules for 
Underwater Technology, 1988 Edition) (reference 3). Design guides included: The Society of 
Naval Architects and Marine Engineers (SNAME) {Submersible Vehicle Systems Design, 1990) 
(reference 4) and the Massachusetts Institute of Technology (MIT) Course 13A Professional 
Summer (13A PS) Submarine Design Trends Course notes (reference 5) was also used. The 13A 
PS is based upon the 1967 version of the SNAME publication Principles of Naval Architecture 
(PNA) with some modifications determined by Harry Jackson (CAPT, USN Ret.). A third 
classification society, the American Petroleum Institute (API) was planned to be used, however 
the immediately available rules from API were not valid for the experimental data used. 
However, API does have rules to cover the types of cylinders examined here. 

Analytic solutions for the three failure modes were gathered from several sources. These 
sources included Hydrostatically Loaded Structures: The Structural Mechanics, Analysis and 
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Design of Powered Submersibles by William Nash (reference 6), Principles of Naval 
Architecture, 1967 edition (reference 1), a David Taylor Model Basin technical paper by J.G. 
Pulos and V.L. Salerno (reference 7), and several journal articles from the Transactions of The 
Royal Institution of Naval Architects byS. Kendrick (reference 8) and C. T. F. Ross (reference 

9 ). 

In Chapter 2 a brief discussion on the terminology of ring-stiffened cylindrical shells is 
produced along with a short derivation of the basic stresses in these cylinders. Chapter 3 
contains the summary of analytic solutions used to predict failure (the solutions are not re¬ 
derived), Chapter 4 has the summary of the design rules and guides from the various sources. 
Chapter 5 describes the test cylinders and then compares the results from experiment to the 
results predicted by the design rules and guides. Chapter 6 summarizes the results with 
recommendations for further study. 



Chapter 2: Overview of Ring-Stiffened Cylindrical Shells 


Ring-stiffened cylinders are the prevalent construction type in the mid-bodies of modem 
submersibles. The stiffeners provide additional strength to the shell that is required for the 
pressure differential between the external hydrostatic pressure and the internal, approximate 
atmospheric pressure. [6] 

2.1 Terminology 

The various classification societies use slightly different terminology for the cylinder 
geometries and properties. In the computerized design rules, the symbols used by each 
classification society are generally used to avoid confusion between the published code and the 
programs. All of the codes require computation of the moment of inertia (I) of the ring stiffener 
and the moment of inertia (I e ) of an effective ring stiffener (the frame and some length of 
attached shell). The formulas for I and I e came from [10]. 

The terminology can be divided into two categories: material properties and geometry. 
All stresses and pressures are in pounds per square inch (psi), lengths are in inches (in), areas are 
in square inches (in 2 ) and moments of inertia are in inches to the fourth power (in 4 ). 

2.1.1 Material Properties 

1) Modulus of Elasticity (E): The slope of the linear region of the stress-strain curve of a given 
material. For all of the calculations a value of 3 X 10 7 psi (a common value for steel) was used. 

2) Poisson’s Ratio (v or p): The ratio of lateral strain to axial strain in a material. For all 
calculations, a value of 0.3 (a general figure for steel) was used. 
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3) Yield Strength (a y or f y or k): An arbitrary value for materials marking the onset of plastic 
deformation of a material. This is usually considered the point of 0.2% permanent strain. This 
value is of primary concern in the axisymmetric yield calculations. This parameter ranged in 
value from 65,500 psi to 157,000 psi. 

2.1.2 Geometry 

2.1.2.1 Cylinder Geometry (see Figure 1) 

1) Cylinder Length (L or Lb or L c ): Overall length of the cylinder between supports. This varies 
by cylinder and is of primary concern for General Instability calculations 

2) Shell Mid-plane Radius (R): Radius from centerline of cylinder to the shell mid-plane. 

3) Shell plate thickness (t or t p or s or h): Thickness of the cylinder shell plating. 

2.1.2.2 Stiffener Geometry (see Figure 1) 

1) Length between Stiffeners (L or Lf or L s ). Distance between centerlines of adjacent stiffeners. 
This distance was assumed to be constant for each cylinder. This dimension is important for 
both Axisymmetric Yield and Lobar Buckling. 

2) Web Height (h w ): Length of the web for ‘T’ stiffeners or the height of the stiffener for 
rectangular stiffeners. 

3) Web Thickness (t w ): Thickness of the web for ‘T’ stiffeners or the thickness of the stiffener for 
rectangular stiffeners. 

4) Flange Breadth (bf): Width of the flange for ‘T’ stiffeners. This value is set to zero for 
rectangular stiffeners. 

5) Flange Depth (d): Thickness of the flange for‘T’ stiffeners. This value is set to zero for 
rectangular stiffeners. 
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6) Stiffener Depth (H): The distance from the shell to the end of the stiffener. This is the 
stiffener height for rectangular stiffeners and the web height plus the flange depth for ‘T’ 
stiffeners. 

7) Faying Length (b): The distance of contact of the stiffener to the shell. It is equal to the flange 
width for ‘T’ and rectangular stiffeners; if T-beam or wide-beam (WF) stiffeners are used, then 
equal to the flange breadth in contact with the shell. 

8) Effective Length of the Shell (L e ): A length of shell to be considered as part of the combined 
stiffener and shell used in the General Instability analyses. This length is usually a fraction of 
the stiffener spacing, but many times will be equal to the stiffener spacing. 

9) Area of the stiffener (A): Cross-sectional area of the stiffener. 

10) Effective Area of the stiffener (A e ff): Cross-sectional area of the combined stiffener and 
effective length of shell or a modified area of the stiffener (A) based on the location of the 
stiffener. 

11) Moment of Inertia (I): Area moment of inertia of the dedicated stiffener, used by some codes 
as part of the calculation for the effective moment of inertia. 

12) Effective Moment of Inertia (I e fi): Area moment of inertia for the combined stiffener and 
shell, used for General Instability calculations. 
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Figure 1: Cylinder and Stiffener Geometry 



2.2 Classification of Stiffeners 

There are several methods to classify stiffeners. Of importance to the analysis of 
cylinders is the location of the stiffener: either internal or external to the cylinder shell. Both 
types are used for cylinder construction; slight modifications to the design rules (concerning 
effective areas and moments of inertia) are required based on the location. External stiffeners of 
equal size to internal stiffeners require more material because of the greater circumferential 
length but offer the advantage of freeing spacing within the cylinder for equipment / living space. 

[ 5 ] 
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Another classification of stiffeners are field (light / non-heavy) and heavy (king frame) 
stiffeners. Heavy stiffeners are substantially larger than the field stiffeners and are used to 
reduce the effective length of the cylinder for general instability concerns. Field stiffeners are of 
uniform size, shape, and spacing. For the current analyses no heavy stiffeners were used. 

2.3 Stresses in Cylinders 

A brief discussion of stresses in cylindrical shells is required to set up the derivation of 
the structural mechanics in support of the analytical solution. As a starting point, a cylinder can 
be considered a thin-walled structure (shell) if the ratio of the mid-plane radius to the wall 
thickness is greater than ten. This assumption allows the determination of the stresses by statics 
alone.[10] All of the cylinders considered are treated as shells. A second assumption in the 
analysis is to consider the hydrostatic pressure as constant across the cylinder. 

Cylindrical shells, exposed to hydrostatic pressure, have two basic stresses imposed by 
the pressure: hoop stress and axial stress. The equations for the stress are: 


1) Hoop (circumferential) Stress: 

^ P R 

h t 

(1) 

2) Axial Stress: 

„ _pR 

*— 

* 2t 

(2) 


Where p is the external pressure, R is the shell mean radius and t is the shell thickness. Figure 2 
shows the derivation of these two equations. 

The addition of ring stiffeners to the base shell complicates the hoop stress analysis by 
introducing non-uniform deformation of the shell in the radial direction. There is also a beam- 
column effect due to the pressure acting in the axial direction. These effects will be addressed in 
Chapter 3. 
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Figure 2: Basic Cylindrical Shell Stresses 



1) Axial Stress 

Axial force on a vertical section of the cylinder must equal the axial stress times the 
circumferential area. 

Force balance: p-tt R 2 = 2-Ji-R t-(T axia| 

Which results in: 

PR 

a axial _ 2 . t 

2) Circumferential (Hoop) Stress 

The force on a transverse section of shell, with width b must equal the hoop stress times the shell 
area. 

Force balance: P-2-R b = 2-b t o hoop 
Which results in: 

PR 

CT hoop •“ t 
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Chapter 3: Analytic Solutions 


Analytic solutions have been proposed and proven for the three major failure modes of 
stiffened cylinders. This chapter briefly describes the theories and equations used in the 
experimental analyses. 

3.1 History of Analyses 

The failure of cylinders exposed to external pressure has been studied for almost 150 
years. The first attempts at understanding cylinder behavior was done by experiment and 
empirical relationships in the 1850’s.[6] The first analytic solution for a non-reinforced cylinder 
was presented by G. H. Bryan in 1888.[9] The first analysis of a reinforced cylinder appeared in 
1913 by R. V. Southwell, followed a year later by a solution to the elastic buckling of a thin shell 
proposed by von Mises.[6] In 1934 Widenburg proposed a solution for asymmetric buckling that 
is independent of the number of lobes of failure. [6] This equation is the one used in the current 
analysis. Solutions for axisymmetric yield were first put forward by von Sanden and Gunther in 
1920.[9] Viterbo presented a modified version of Sanden and Gunther’s solution in 1930. [9] 
Pulos and Salerno presented a solution that included the Sanden and Gunther solution, the 
Viterbo modification and a term to account for the bending stress in the cylinder caused by the 
axial pressure.[9] The Pulos and Salerno solution is used in this thesis. For elastic general 
instability, the first reported analysis was presented by Tokugawa in 1929. In 1954 A. R. Bryant 
developed a similar equation using a different methodology. [6] 

Analytical work from the 1950’s onward has focused on obtaining solutions for different 
boundary conditions and more fully reconciling the analytic predictions with experimental 
results and to more fully understand the effects of initial imperfections in the cylinder’s material 
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and geometric properties. With the advent of power digital computers and the use of finite 
element analysis there has been great strides made in understanding the failure of cylinders. 

3.2 Current Theory 

For this thesis, analytic solutions were collected from several sources for the modes of 
failure. A comprehensive theoretical solution that addresses all modes of failure is not presented. 
Reference 6 provides a good summation of the currently used analytic solutions. 

A first indicator of the failure mode of a cylinder is found by plotting the cylinder’s 
slenderness ratio (X) against the ratio of the shell buckling pressure(p c ) to the hoop pressure at 
yield(py) (v|/).[l ] The equations for these factors are given below. 


i 



Figure 3 shows the plot of \\i versus X. If the slenderness ratio is less than approximately 1.14 the 
cylinder should fail by axisymmetric yield; if l is greater than 1.14 then the cylinder should fail 
by lobar buckling. If the shell and stiffeners are not of sufficient size, the cylinder may fail by 
general instability at a pressure less than that found in Figure 3.[ 1 ] By using two assumptions 
(the material is steel with t> = 0.3 and L s /D » t p /D) it can be shown that: [1] 

1.30 

< 5 > 

This is the buckling part of the curve in Figure 3. 
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Figure 3: Failure Pressure Ratio versus Slenderness Ratio 



A 


The treatment of boundary / end conditions of the cylinder is a vitally important factor in 
the analytic solutions. The literature is full of discussion on what types of end condition to use, 
with the choices ranging from full clamped ends to simply supported ends. The extreme cases 
are hard to create in reality and therefore the experimental results tend to fall between the ranges 
of predictions. Experiments have shown that partially clamped cylinders provide significantly 
higher failure pressures than that predicted by mathematical models utilizing simply supported 
ends.[8] For this thesis, no discrete boundary conditions were required to be stated for input into 
the equations. 

3.2.1 Axisymmetric Yield 

As mentioned in section 3.1 axisymmetric yield has been studied since the 1920’s. The 
solution summarized here was put forward by Pulos and Salerno in 1961. It is based on the 
previous works of van Sunden and Gunther and Viterbo and includes a previously not included 
“beam - column” effect due to the hydrostatic pressure acting in the axial direction of the 
cylinder.[7] The governing differential equation is: 
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( 6 ) 



r \ 




D 

d 4 

vv 

+ 2±. 

d 2 Et 

-H) 


A 4 

^ dx J 

2 

- \y -f.-\y 

dx 2 R 2 


Where: 

Et'* 

D := —-r D is the flexural rigidity of the shell 

12-ii - v 2 J 

The term represents the beam - column effect. It makes the solution to equation (6) a non¬ 
linear function of pressure and was the term neglected in the earlier analyses of axisymmetric 
yield. Figure 4 shows the coordinate system used in reference 7 to derive the governing 
equation, x, <p, and r are the axial, circumferential and axial coordinates respectfully with u, v, 
and w being the corresponding displacements. [7] 


Figure 4: Coordinate System for a Cylindrical Shell 
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Following typical practice in the solution of non-homogeneous differential equations the 
general solution of the governing equation can be written as the sum of the solution of the 
homogeneous equation and a particular solution. [6] The solution to the homogeneous equation 
produces four roots (X1-X4). By analysis, placement of the origin of the coordinate system to take 
advantage of symmetry, and trigonometric identities, the general solution can be given as: 

w — B cosh + E cosh T, x — — (1 — —) (7) 

Et 2 

where B and F are arbitrary constants which can be found by applying boundary conditions to 
the equation.[7]. After further mathematical substitutions several dimensionless parameters were 
introduced into the solution to allow easier solving of the problem. Four of these parameters (Fj- 
F4) were transcendental functions based on the geometry of the cylinder. These functions were 
originally graphed to allow for a relatively quick solution to be found for a given cylinder. 

Finally an equation for the failure pressure of a given cylinder was determined. This equation 
along with the dimensionless parameters is given below. 
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1-Ua 


P + 0 - P)-F, 


Y:= ^E 3 1 ~ 


0 , := - ypT~y 


*i2 := -V^ r r 


0 := \/ 3 ll - |i 


2 L 


2 2 

cosh (r| j-0^ -cos ( r >2 0 ) 
cosher] |- 0 ^*sinh^T| j- 0 j cos^^-sin (0 2' 6 ) 

01 02 

-sh^r| j-0ysin^ri2-0j sinh^rj j-0jcos^r^*©) 

^2 _Hj_ 

sh ^rj j*0 j-sinh^rj j-0^ cos^2-0^-sin (0 2 e) 

01 0 2 

cos (0 2' 9 ) •sin^02'0j cosh^0 j-0^-sinh|0 j-0^ 

3 _ 02 _ 01 _ 

2 cosh^0 j-0^ sinh^0 j-Oj cos^-O^-sin^jO) 

01 02 

cosh^0 |-0j-sin^02-0) sinh^i j-Gj-cos^-O) 

3 _ 02 _ 0 j_ 

_ ^2 cosh^0 |-0j sinh^0 j-0^ cos^-Gj-sin (0 2’ 0 ) 


An iterative process is required for the general case where the parameter y is not zero. The 
process begins with assuming that y is zero, finding the failure pressure then recalculating y and 
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solving the equations for the failure pressure again. Usually only two to three iterations are 
needed for satisfactory convergence of the failure pressure. [7] 

3.2.2 Asymmetric Buckling (Lobar Buckling) 

Asymmetric buckling is characterized by circumferential lobes between ring stiffeners. 
As noted above this mode of failure will occur when the cylinder’s slenderness ratio is relatively 
high. This can be further characterized by a relatively thin shell thickness and widely spaced 
stiffeners.[l] R. von Mises first proposed a solution to the buckling of un-stiffened cylinders 
under hydrostatic pressure in 1929. He assumed sinusoidal displacements in the axial and 
circumferential directions to enable solving of a set of linearized partial differential equations 
that represented the elastic action of the shell. [6] He eventually obtained the following equation 
for the buckling pressure: 



The buckling pressure is dependent upon the number of circumferential lobes (n), which must be 
an integer value. This fact requires an iterative process of varying n until the lowest pressure is 
determined. 

In 1933 Widenburg solved the above equation in way that was independent of n. From 
test data the buckling pressures from the von Mises equation and the Widenburg approximation 
differ by no more than 3.5%. [6]. Further investigation into the bucking of stiffened cylinders 
determined that the Widenburg equation worked very well by replacing the length of the cylinder 
by the length between stiffeners. Therefore the Widenburg equation (equation 10), shown 
below, is the equation used in the current analyses. 




3.2.3 General Instability 

General instability consists of the yielding of both the cylindrical shell and ring stiffeners. 
A cylinder may be susceptible when the stiffeners are undersized when compared to the shell 
thickness and the cylinder is relatively long.[6] General instability may occur in either the elastic 
or inelastic stress region of a material. Elastic general instability is the mode covered by the 
available literature and is addressed here. Inelastic general instability has been studied mainly by 
government laboratories and most of the knowledge of this failure mode is classified material 
and therefore unavailable to the present author. [11] 

The first analysis of general instability was presented by T. Tokugawa in 1929. [6] His 
methodology considered the failure as a combination of the failure of the ring stiffeners and shell 
buckling with each taking place separately. [9] In the 1940’s S. Kendrick used a strain energy 
methodology to determine the failure pressure. In 1953 A. R. Bryant used a simpler strain 
energy method and determined an equation similar to both Kendrick’s and Tokugawa’s. [1, 9] 
The Bryant equation is therefore used for this analysis and is shown below as equation 10. 
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Lb 


The first term can be considered the failure of the shell and the second term can be 
considered the failure of a combined stiffener and an effective length of shell. [6] This effective 
length has had much discussion in the literature over the decades. Bryant assumed the length to 
be the spacing between the stiffeners, but others have proposed various corrections based on the 
cylinder’s geometric and material properties.[5] For this thesis, the effective length used was the 
stiffener spacing. 

Similar to the von Mises buckling pressure determination, the number of circumferential 
lobes (n) must be found that minimizes the failure pressure. The number of circumferential lobes 
is usually between 2 and 4. [5] The factor X is the number of longitudinal lobes in the cylinder. 
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Chapter 4: The Design Rules and Analysis Tools 


The two classification societies’ design rules that were utilized were the American 
Bureau of Shipping and the Germanischer Lloyd rules. These were chosen for their availability, 
their different levels of simplification of equations, their coverage of the specific geometries of 
the experimental cylinders, and their inclusion of all three failure modes. The Society of Naval 
Architects and Marine Engineers Submersible Vehicle Systems Design and Principles of Naval 
Architecture (PNA) were chosen for analysis as SNAME is the primary design society in the 
United States and has comprehensive guidelines for cylinders. The MIT 13A Professional 
Summer Submarine Design Trends notes were used because of the author’s familiarity and the 
complete analysis of the failure modes. 

4.1 American Bureau of Shipping (ABS) 

The ABS design rules, as delineated in the Rules for Building and Classing Underwater 
Vehicles, Systems and Hyperbaric Facilities (1990), give a brief and conservative approach for 
determining the critical / collapse pressures for each failure mode. The ABS design rules do not 
explicitly name the failure modes, but there are distinct equations for the three modes. 

4.1.1 Axisymmetric Yield 

This mode is designated the yield pressure at midbay and midplane of a cylinder. [2] The 
formula accounts for the major parts of the analytic solution, but uses single value functions for 
the shell parameters and does not explain the functions of each part of the equation. The rules do 
account for the difference between internal and external frames by squaring the mid-plane radius 
to stiffener radius ratio for external stiffeners. The base equation follows: 
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( 12 ) 


t 



Where 


F := 



A + t w -t + 


2NtL 

6 



Internal Stiffeners or 



2 

External Stiffeners 



N — cosh ( 2 'Q) ~ cos(2 Q) 
sinh(2-Q) + sin(2Q) 

G = : (sinh(Q) cos(Q) + cosh(Q) sin(Q)) 
sinh(2Q) + sin(2Q) 

H . sinh(2Q) - sin(2 Q) 
sinh(2-Q) + sin(2Q) 


4.1.2 Asymmetric Buckling 

This is called the Von Mises buckling pressure for a cylinder. It is the Widenburg 
approximation that is used by most of the classification societies. 
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The ABS code then has a range of allowable pressures depending upon the ratio of P m to 
P y . Below is the logic for the maximum allowable working pressure for inter-stiffener strength 
(P c ). Further safety factors are then applied to lower the allowable pressure in practice. 


P P 

m m 

— if — <1 
2 P 

y 

( p \ p 

y m 

P • 1-— if 1 < — <3 

y 2-p p 


5 .. m 

-P if — >3 

6 y P y 


4.1.3 General Instability 

For the elastic general instability, the ABS code uses the Bryant equation with an 
effective length not equal to one frame space. The effective length is given below. 


L := mini 1 5 The equation for the failure pressure is broken into three parts, but the 

e I 075L S 


total is equivalent to the Bryant formula. 


Et El e A 2 (n) 

Pnl< n ) :=—•A 1 (n) +--- 

R r 3 l 


Where: 
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4.2 Germanischer Lloyd 

The design rules for Germanischer Lloyd are from Rules for Underwater Technology 
(1988 edition). The rules address all three failure modes very thoroughly and flow charts are 
provided to aide in programming the code for computer use. These rules also address out-of- 
roundness up to a nominal value of 0.5% of the mid-plane shell radius by determining a reduced 
allowable pressure.[3] 


4.2.1 Axisymmetric Yield 

The code for axisymmetric yield resembles the Pulos and Salerno methodology with 
some additions to account for the transition of the material into the plastic range during yield. 
The methodology consists of first guessing a pressure lower than the failure pressure (the test 
cases used 1 psi as the starting point) then iterating through a series of equations to determine the 
shell stresses. If the determined shell stress becomes greater than 0.8 of the yield strength of the 
material, then a Secant Modulus of Elasticity (E s ), a Tangent Modulus of Elasticity (E t ) and a 
plastic range Poisson s Ratio (Up) is calculated and used to determine the other calculation 
factors. The determination of the integer m in equation (14) is described in the code below. 







h 3 h 4 


C3 !— 1 + 


h 2 h 3 h 4 


C 4 is not used for cylinders without heavy stiffeners. 


C5a-Lj 


Past • 


2-s 2 E 

R 2 J 3(1 v 2 ) 


critical pressure, elastic, calculation factor 


p = 1 psi (arbitrary low pressure) 

This is the actual starting point of the iteration. A pressure (p) is guessed; the code then 
calculates the failure pressure. If the failure pressure is greater than the guessed pressure, the 
applied pressure is increased by a set amount. This continues until the two pressures, p and pm 
are within a chosen delta of each other. This is essentially the method of Pulos and Salerno with 
a simpler method for iterating on the presumed pressure. 


C 6 := - VT -G 


C7— i/T+G 


C 8 C 5' C 6 


C 9 : “ C 5 C 7 


i-n-A. 

2) s-L, 


A b ( b) 
h h 1-Fi 

SL 1 L 1 L 1 


32 



c ll : “ 


0.91 


1 -v 


Fi := —■ 


2 2 
cosh^Cg) -cos^Cp) 


^ cosh(Cg)-sinh(Cg) cos(Cp)-sin^Cp) 


cosh(Cg)-sin^Cp) sinh ^Cg)-cos (Cp) 


F 2 ^ 


cosh (Cg)'Sinh^Cg) cos •sin N 


cos (Cp) ■ sin^Cp) cosh (Cg) - sinh (Cg) 


F 3 := 


2 cosh (Cg)-sinh (Cg) cos N •sin^Cp^ 


cosh ^Cg^ • sin^Cp) sinh (Cg)* cos ^Cp) 


F 4 := 


2 cosh (Cg)-sinh(Cg) cos(Cp)-sin(Cp) 

1 ^ - _i_- 


The functions FI through F4 are equivalent to those presented in the analytic solution. 


a := 




s 2 -R 2 


shape factor 


The following set of equations determines the stresses at the mid-bay of the shell. 


a o := 


-p-R 


°x : -V - + c 10' c 11' F 4 


(7 <|) : - a o'( 1 “ C 10' F 2 + v C 10' C ir F 4) 


C xl : -V 2 _C 10' C lf F 4 


°(j)l “ C 10 F 2 " v C 10' C l r F 4) 
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If the calculated stress is > 0.8k then the following equations are used to determine the 
material properties in the elastic - plastic region. The code is somewhat circular, the assumption 
was made to determine the strain (ei) from the previously calculated stress then to find a new 
stress level by the equation given in the code. 
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4.2.2 Asymmetric Buckling 


For asymmetric bucking the Germanischer Lloyd code uses a modified version of 
equation 9. Therefore the dependence upon the number of lobes is not removed and should give 
a better estimate of the buckling pressure. The equation is given below. 


PniW : 


Es-P nl (n) 


Where: 


Pnl(") ; =- 


2 y 2 2 f 2 . . 2 

n 2 , s -[n - 1 + A.j 

- 2 -2.R 2 (,-v 2 ) 


2 , ..2 
n — 1 + .5‘X| 


4.2.3 General Instability 

For elastic general instability, the code uses the Bryant equation with the effective 
moment of inertia based upon the effective length of shell as defined: 

L etest := b + V 2 ' R s 

L e := L etest if L etest < Lj where Li is the frame spacing 
Lj otherwise 

The base equation for general instability also uses a modified radius (Ro) which is the radius to 
the centroid of the combined stiffener and effective length of shell. For this thesis no “heavy” or 
“king frame” stiffeners were used, however the code allows for these stiffeners with several 
more equations. 


PnlaW := P 0 ( n ) + Pl( n ) 
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Where: 



4.3 Society of Naval Architects and Marine Engineers (SNAME) 

The SNAME code was taken from Submersible Vehicle Systems Design (1990). SNAME 
is not a classification society as ABS and Germanischer Lloyd, but it does provide many of the 
professional and technical resource material for naval architects in the United States. The 
current codes are revised versions of those found in the SNAME publication Principles of Naval 
Architecture (1967 edition), reference 1. The code addresses all three modes of failure with the 
least complex set of equations of those under current study. 

4.3.1 Axisymmetric Yield 

For axisymmetric yield, the methodology simplifies the analytic solution by using the 
older theory of von Sanden and Gunther instead of the solution put forward by Pulos and Salerno 
in reference 7. The methodology thus neglects the beam-column effect of the axial pressure. 

This is justified by designing the shell to yield vice buckle.[1] 
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4.3.2 Asymmetric Buckling 

For Lobar buckling, the code uses the Widenburg approximation equation that is 
independent of the number of failure lobes. This is the same as used by ABS. 
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4.3.3 General Instability 


For general instability, the code uses the unmodified Bryant equation; therefore the 
length of shell used in the moment of inertia calculation is one stiffener spacing. The code 
suggests that the critical value of n is between 2 and 4. 


Pcrl( n ) := 


E-t 

R 


m 


(n 2 - i) e-L 


f 2 ^ 

2 t m 
n - 1 + — 

2 


( 2 2 ) 
\n + m / 


R Lr 


( 20 ) 


Where: 


m := 


7C- R 

h 


4.4 MIT 13A Professional Summer Submarine Design Trends 

The Professional Summer notes are a compilation of design theories from Harry Jackson 
(CAPT, USN, Ret.) used to instruct the MIT course 13A students in the basics of submarine 
design. The actual purpose of the code is to evaluate a proposed submarine design given depth 
and material criteria; some modifications were required to provide the failure pressures. The 
failure equations for lobar buckling and general instability are generally from reference 1, 
however some small modifications have been made. A complete stress calculation is included 
for axisymmetric yield. 

4.4.1 Axisymmetric Yield 

The majority of the modifications to the original code were made to allow calculation of 
a failure pressure for this mode. The beginning code provided a required diving depth and 
therefore a pressure to withstand. It then performed a von Mises stress analysis on the structure, 
resulting with a shell stress. This shell stress (after modification by a safety factor) was then 
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compared to the yield strength of the material for acceptability. The code was modified to iterate 
on the pressure, with the failure pressure designated as the pressure at which the calculated stress 
was equal to the yield stress (without any modifications by safety factors). There was no distinct 
equation for the failure pressure. The following is the code used to determine osy, which is the 
shell stress, 
hv'tp 

B :=-— 

A + tvv-t P Area ratio 


l 



cosh ( 9 ) - cos ( 9 ) 
sinh(9) + sin(9) 


Slenderness Parameter 


Deflection coefficient 


P:= 


2-N 

A + t w -tp 



Frame flexibility parameter 



Frame deflection parameter 


Bending effect (mem) 


H E :=-2- 


3 

2 

1 -v 



sinh(9) - sin(9) 
sinh(9) + sin(0) 


The following begins the von Mises stress analysis. 


Bending effect (bend) 


Bending efffect near frame 
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4.4.2 Asymmetric Buckling 

The equation for lobar bucking is the same as that used by SNAME and ABS and the 
analytic solution; it is the Widenburg approximation: 
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For general instability, the Bryant equation is used with the effective length of shell given 
by the following method: 

Lgff :=L C F] + t w 
Where 

Lg := Lf - t w is defined as the clear length between the mid-bay and the stiffener 

2 2 

4 cosh(nj-0) - cos(n2*0) 

P j ♦ — —. - 

0 cosh(n]'0)sinh(nj-0) cos(n2'0)*sin (n 2 e) 
n l n 2 

nj := 0.5 t/ 1 - y n 2 := 0.5-yj 1 + y 



F) is the first transcendental equation determined by Pulos and Salerno. The original 
methodology required the pressure as an input to first determine the 0 function (it is the same as 
for axisymmetric yield) and then Fi was found. To provide for a simpler method, a series of 
pressures that encompassed the predicted range of failure were input into the above equations 
and an average value of Fi was chosen from the resultant values. The Fi values changed by 
approximately + 0.05 for each test cylinder. The effective length was then used to calculate the 
moment of inertia for the combined shell and stiffener. The number of lobes was set to be 
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between 2 and 4; the failure pressure was then found without other modifications to the Bryant 
formula. 


p cGI : = 

l 


R 





E-Icff 


R 3 Lf 


( 22 ) 


Where: 



Chapter 5: Experimental Results 


The various Classification Societies’ design rules were tested against data collected from 
experiments conducted by the U. S. Navy in support of submarine design. Each society’s design 
rules were used to determine the failure pressure and failure mode of each of the test cylinders. 
The resulting predictions were then compared to the experimental results. Of primary interest 
was agreement between the design rule prediction and experiment on the mode of failure, 
followed by the closeness of the predicted failure pressure to the actual collapse pressure. 

5.1 NAVSEA Test Cylinders 

The test data was provided by the Naval Sea Systems Command submarine structures 
unit (SEA 05P2). Four test cylinders were selected that covered the range of examined failure 
modes. The cylinder diameter to thickness ratios (D/t) fell between 112 to 198, modeling typical 
submarine D/t ratios. Two of the cylinders had internal stiffeners while the other two cylinders 
had external stiffeners. All four test cylinders had built-up end stiffeners with a combination of 
slightly different spacing and / or larger stiffener dimensions than the uniform field stiffeners. 
These end stiffeners were designed to prevent shell yielding in the end bays due to increased 
stress levels associated with the boundary conditions. It was estimated that without the end 
stiffeners a 4-5% reduction in axisymmetric yielding pressure could occur. None of the end 
stiffeners met the classification societies’ requirement for a “deep” stiffener. None of the design 
rule codes allowed for these variable stiffeners, therefore the non-uniformities were disregarded 
and the end stiffeners treated as field stiffeners. The four cylinders are described below. 
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5.1.1 Cylinder l.d 


Cylinder 1 .d was a machined cylindrical shell with rectangular external ring stiffeners. 
The material was high strength steel with a yield strength of 80,000 psi. Figure 5 shows the 
schematic of the cylinder. The boundary conditions consisted of one end being fully fixed; the 
other end had all freedoms fixed except for axial displacement (these conditions conflict with the 
design rules assumption of completely clamped ends). External hydrostatic pressure was applied 
including axial line load to simulate load on the end plate. The experiment tested the ability of 
the design rules to predict elastic shell bucking (Lobar buckling). The experimentally 
determined collapse pressure was 633 psi with failure by asymmetric (Lobar) buckling. 
Appendix A has the analysis of predicted failure for this cylinder. 

Figure 5: Test Cylinder l.d Schematic 
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5.1.2 Cylinder l.f 


Cylinder 1 .f was a cylindrical shell with internal tee stiffeners of welded construction. 

The material was high strength steel with a yield strength of 98,500 psi. The boundary 
conditions consisted of 4.0 inch steel plates attached with full fixity to the end of the adaptor ring 
on the model. External hydrostatic pressure was applied. This test cylinder was used to predict 
failure by elastic general instability. There was no experimental elastic collapse pressure; 
therefore the critical pressure was calculated by two separate, reliable analysis programs with the 
results being 4858 psi (with 3 waves) and 4953 psi (with 3 waves). The test cylinder actually 
failed by inelastic general instability at a pressure of 2200 psi. Figure 6 shows the cylinder 
dimensions and Appendix B has the analysis of predicted failure for this cylinder. 


Figure 6: Test Cylinder l.f Schematic 












5.1.3 Cylinder 2.a 

Cylinder 2.a was a machined cylindrical shell with external tee stiffeners. The material 
was high strength steel with a yield strength 65,500 psi. Figure 7 shows the schematic of the 
cylinder. The boundary conditions consisted of end closures made of 3.0 inch steel plates 
attached to the idealized adaptor ring with full fixity. External uniform hydrostatic pressure was 
applied to the model. This cylinder was used by the Navy to predict end bay failure (shell 
collapse influenced by end bay design). This is a specific example of axisymmetric buckling and 
was used as the axisymmetric model for the classification society rules. The experimental 
collapse pressure was found to be 921 psi by axisymmetric collapse in the second bay from the 
adaptor ring. Appendix C has the analysis of predicted failure for this cylinder. 

Figure 7: Test Cylinder 2.a Schematic 
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5.1.4 Cylinder 2.c 


Cylinder 2.c was a fabricated cylinder with internal ring stiffeners. The base material 
was high strength steel with a yield strength of 157,000 psi. Figure 8 shows the schematic of the 
cylinder. The shell was cold rolled and fabricated with a deliberate out-of-roundness 
imperfection. The frames were built-up. The frame web material was base metal, and the frame 
flanges were cold rolled. The boundary conditions consisted of one end being fully fixed with 
the other end having all freedoms except axial displacement. External uniform hydrostatic 
pressure with an axial end load to simulate end plate loading was applied. This test cylinder was 
used by NAVSEA to predict the inelastic general instability failure mode and to model out-of- 
roundness imperfections. In the current comparison the out-of-roundness was disregarded. The 
collapse pressure was experimentally found to be 3640 psi in 2 circumferential waves in an 
inelastic general instability mode. Appendix D has the analysis of predicted failure for this 
cylinder. 

Figure 8: Test Cylinder 2.c Schematic 








5.2 Calculation to Experiment Comparison 


Table 1 shows the comparison of the classification societies’ design rules to the 
experimental results. The table displays the calculated failure pressure and failure mode for each 
cylinder and the percent difference from the experimental failure pressure (if the experimental 
and calculated failure modes are the same). 


Table 1: Comparison of Design Rule Calculations to Experimental Results 



1 .d 



2.a 

2.c 1 


Pressure 

Mode 



Pressure 

Mode 

Pressure 

Mode 

NAVSEA 

EPERIMENT 

633 

■ 

2200 

iGI 

921 

AX 

3640 

iGI 

Analytic 

Solution 

605 

-4.6% 

■ 

2141 

AX 

876 

-5.1% 

AX 

4080 

AX 

ABS 

605 

-4.6% 

■ 

2039 

AX 

844 

-8.8% 

AX 

4211 

AX 

PNA 

605 

-4.6% 

■ 

1928 

AX 

815 

-12.1% 

AX 

3864 

AX 

Germanischer 

Lloyd 

606 

“4.5% 

■ 


AX 

1030 

12.4% 

AX 

4567 

AX 

13A Professional 

Summer 

605 

“4.6% 




819 

-11.6% 

AX 

3712 

AX 


Key: 

L 

Asymmetric (Lobar) Buckling 


AX 

Axisymmetric Yielding 


GI 

General Instability 


e/i 

elastic / inelastic 


There is excellent agreement between the experimental data and the calculations for 
cylinder 1 .d. The Lobar buckling failure was expected as the slenderness ratio was 201, several 
magnitudes greater than the breakpoint of 1.14 between asymmetric and axisymmetric failure. 
The agreement between the design rules and the analytic solution calculations was expected as 
they all used the same formula to determine the lobar buckling pressure. The calculated failure 
pressure was 4% below the experimental pressure. The higher experimental critical pressure was 
attributable to the test cylinder being more fully clamped than theorized in the design rules. 
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Cylinder 2.a also generated agreement between the experimental data and the design rule 
calculations. The predicted failure mode of axisymmetric yield was confirmed by the 
experiment. However there were some significant differences between the predicted collapse 
pressure and the experimental failure pressure. All of the design rules, except for the 
Germanischer Lloyd, predicted failure at a pressure lower than that experimentally found. The 
Germanischer Lloyd calculation over estimated the failure pressure and was the furthest from the 
experimental pressure with a pressure 11.8% over the experimental pressure. This failure mode 
is the most dependent upon the yield strength of the material therefore a small variance between 
the given yield strength and the actual yield strength of the test specimen may have contributed 
to the differing pressures predicted by the design rules. All of the design rules use a simplified 
version of the methodology presented by Pulos and Salerno [7]. The Germanischer Lloyd and 
the modified Professional Summer calculations, perform an iterative operation to find the 
pressure, very close to the analytic methodology. The ABS and PNA calculations use single 
value equations substituted for the transcendental functions of the analytic solution. 

Cylinder 1 .f failed at 2200 psi experimentally in an inelastic general instability mode. 

The experiment was done to test the ability of NAVSEA’s computer codes to predict elastic 
general instability. A predicted failure pressure of 4858 psi was determined for the elastic 
general instability mode. The design rule codes estimated the general instability failure pressure 
very well. Table 2 compares the design rules and analytic solution general instability pressures 
to the experimental failure pressure. The agreement between the design rules was expected as 
they all use the Bryant equation (equation 10) to determine the failure pressure. The only 
differences come from the variations in the effective length of the shell for the combined shell 
and stiffener calculations and small variations in the radius used in the equation. 
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Of further interest to cylinder 1 .f was that the design rules predicted failure by 
axisymmetric yielding at an average of 2206 psi.. As the cylinder actually failed by inelastic 
general instability the design code predictions indicate that the two modes of failure are very 
close together. This cylinder was very close to failure in multiple modes at approximately the 
same pressure (i.e. an identical cylinder of the same dimensions and material may have failed by 
axisymmetric yielding vice the general instability depending on the eccentricity of the cylinder 
and other small defects). This multiple failure mode condition must be guarded against in real 
designs, usually by applying different safety factors to the various modes. [1] 

Table 2: Cylinderl.f Elastic General Instability Failure Pressures 


DESIGN RULE 

FAILURE PRESSURE 
(PSI) 

% FROM 
EXPERIMENT 

LOBES 

Experiment 

4858 

— 

3 

Analytic Solution 

4496 

- 7.5 

3 

ABS 

4496 

- 7.5 

3 

PNA 

4496 

- 7.5 

3 

German!scher Lloyd 

4651 

- 4.3 

3 

13A PS 

4460 

- 8.2 

3 


Cylinder 2.c failed at 3640 psi in an inelastic general instability mode with 2 
circumferential waves. However the design rule calculations all predicted failure by 
axisymmetric yielding at an average pressure of 4086 psi. For this cylinder the design rule codes 
were not close in predicting the elastic general instability failure pressure. This large 
overestimation of the failure pressure can be attributed to the assumption of perfect circularity in 
the design codes, whereas the experimental model had a two-wave sinusoidal imperfection of 
maximum height of + 0.105 inches. This deliberate out-of-roundness would significantly reduce 
the resistance to buckling. Table 3 compares the design rules / analytic solution to the 
experimental failure pressure. 
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From Table 3, the axisymmetric failure pressures were higher than the actual failure 
pressure but all were within 16% of the experimental result. This closeness between the failure 
modes resembles the results for cylinder l.f. Further study may be warranted to explore 
connections between axisymmetric yield and inelastic general instability. 


Table 3: Cylinder 2.c Elastic General Instability Failure Pressures 


DESIGN RULE 

FAILURE PRESSURE 
(PSI) 

% FROM 

EXPERIMENT 

LOBES 

Experiment 

3640 

— 

2 

Analytic Solution 

8642 

137.4 

2 

ABS 

8642 

137.4 

2 

PNA 

8642 

137.4 

2 

Germanischer Lloyd 

9702 

166.5 

2 

13A PS 

8536 

134.5 

2 i 
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Chapter 6: Conclusions 

The classification society design rules studied are important tools for engineers and naval 
architects designing and studying cylindrical structures subjected to external hydrostatic 
pressure. The engineer must have confidence that the design code used will provide acceptable 
(safe) calculations for his or her structure. This confidence can be assured by comparison of the 
calculated failure pressure and mode to that found from experiments. This thesis attempted to 
provide that comparison for several of the most used design rules along with a comparison of the 
analytical solution upon which most the design rules are at least partly based. 

6.1 Comparative Analysis Review 

As discussed in Chapter 5, the design rules had mixed results in correctly predicting the 
failure mode and failure pressure of the test cylinders. In general all of the design rules and the 
analytic solution were in agreement for the specific geometries. For the test cases of failure by 
axisymmetric yielding and lobar buckling, the calculated pressures were accurate when 
compared against the experimental results. However for the two cylinders that experimentally 
failed by general instability, all of the design codes predicted failure by axisymmetric yielding 
vice the general instability. The design rules only account for elastic general instability which 
will occur at a higher pressure than the inelastic general instability. For cylinder 1 .f the design 
codes calculation of the elastic general instability mode failure pressure was very accurate, but 
for cylinder 2.c the calculated failure pressures were all greater than 100% over the experimental 
pressure. The majority of the pressure differentials for 2.c can be attributed to the built-in non¬ 
circularity in the test cylinder. 
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6.2 Agreements and Differences 


In this thesis there were thirty failure pressures calculated (five failures for each of four 
test cylinders and five general instability pressures each for the two cylinders that failed by 
general instability) with design rule / analytic solution agreement on 33% of the pressures. The 
other twenty pressure calculations varied from over 30% to less than 1% different. 

The pressure calculations that were in agreement were generated from the use of the same 
equation and same dimensions for most of the asymmetric buckling and elastic general 
instability calculations. For the asymmetric buckling predictions, ABS, PNA, and the MIT 13A 
Professional Summer calculations use the same equation as the analytic solution. This equation 
was independent of the number of circumferential lobes in the failed part of the shell. For the 
general instability cases, the ABS and PNA calculations agreed with the analytic solution. These 
three codes used the same effective length (L c = frame spacing) and the assumption that the 
radius of the combined shell and stiffener was the mid-plane radius of the shell. 

The difference between calculations can be accounted for individually: The 
Germanischer Lloyd code for asymmetric buckling used an iterative process that was dependent 
of the number of lobes in the failed cylinder; however the failure pressure was still within a tenth 
of a percent of the other calculations. The predictions for elastic general instability were also 
generally close. The Germanischer Lloyd and MIT 13A Professional Summer codes used 
different assumptions for L e than the other codes, which gave predicted failure pressures slightly 
different than the base general instability equation (equation 10). The largest variations between 
the codes were generated for the axisymmetric yield mode. All of the codes used simplified 
variations the analytic solution developed by Pulos and Salerno. These variations in the shell 
yielding pressure come about due to the different simplifications made in the codes. 
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6.3 Applications of the Models 

The various design rules studied are promulgated to ensure safe design of cylinders for 
use under external pressure conditions. This safe operation requires an almost absolute certainty 
that the design will not fail under the worst anticipated condition (many theories and practices 
exist on risk based design and the use of safety factors). The comparison of the design rules to 
an analytic solution, analysis tools (MIT 13A Professional Summer and SNAME), and to 
experimental results allows a designer to have a good idea of how a particular cylinder design 
would be evaluated by each entity. This comparison ability would be useful in judging the initial 
feasibility of a design and also would be useful in applications where a design would be subject 
to more than one classification society. 

The comparisons in this thesis should in no way be used as a detailed design tool for the 
subject cylinders. After an initial design is compared and judged to be adequate, much more 
rigorous analyses must be used to ensure a safe design. These advanced analyses should include 
finite element methodology and other tools that can look at local stresses instead of generalized 
stresses in the shell. These higher order analysis tools can account for material differences, 
geometric eccentricities (out-of-roundness and other along the shell), varied spacing and sizes of 
stiffeners, and actual construction factors such as heat effected zones around welds, and 
bulkhead effects. Another important area that is addressed by other parts of the design rules but 
not studied here is that of shell penetrations. These discontinuities in the shell are very 
susceptible to stress concentration and must be reinforced to prevent failure at pressures lower 
than that predicted for a continuous shell. 
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6.4 Further Areas of Study 

There are several areas that require further research to completely understand the result of 
the above analysis. The area most evident in need of more study is the failure of the cylinders by 
inelastic general instability. While some work has been accomplished on this phenomenon, most 
of the work is either classified by government entities or is empirical data. As part of the 
research into the inelastic general instability should be a detailed analysis of interactions / 
relationships with axisymmetric yield. As found from the experimental cylinders that failed by 
inelastic general instability, the design codes and analytic solution predicted failure by 
axisymmetric yield of the shell at pressures close to the actual failure pressures. These 
comparisons suggest that at least the two modes are very close together and may have some 
interaction. 

There are many other classification societies that produce design rules for stiffened 
cylinders and other geometries. These societies include the American Petroleum Institute (API), 
NORSOK and Det Norske Veritas (Norway), Lloyd’s Register (United Kingdom), Registro 
Italiano Navale Group (RINA) (Italy) and several others. These additional design rules could be 
compared against the existing test cylinders as well as all of the design rules should be compared 
against more experimental data. 

Finally, the base methodology can be expanded to include other geometries such as 
spheres, hemi-spheres, cones and toroids. The classification societies address most of these 
forms which can be found as the end closures for most submersibles and other pressure vessels. 
These more complex geometries present more challenging analytic solutions and are in general 
harder to manufacture. Experimental data is scarcer; therefore finite element models may be 
needed to generate comparisons. 


56 


References 


[1] Society of Naval Architects and Marine Engineers. 1967. Principles of Naval 
Architecture. New York: The Society of Naval Architects and Marine Engineers. 


[2] American Bureau of Shipping. 1990. Rules for Building and Classing Underwater 
Vehicles, Systems and Hyperbaric Facilities. New York: American Bureau of Shipping 

[3] Germanischer Lloyd. 1988. Rules for Underwater Technology. Hamburg: Germanischer 
Lloyd 

[4] Society of Naval Architects and Marine Engineers. 1990. Submersible Vehicle Systems 
Design. Jersey City: The Society of Naval Architects and Marine Engineers. 

[5] Massachusetts Institute of Technology Program 13A Professional Summer. 2001. 
Submarine Design Trends. Unpublished 

[6] Nash, William A. 1995. Hydrostatically Loaded Structures: The Structural Mechanics, 
Analysis and Design of Powered Submersibles. Tarrytown: Pergamon Press. 

[7] Pulos, J.G., and Salerno, V.L. 1961. Axisymmetric Elastic Deformations and Stresses in a 
Ring-Stiffened, Perfectly Circular Cylindrical Shell Under External Hydrostatic Pressure. David 
Taylor Model Basin Report 1497. 


[8] Kendrick, S. 1965. The Buckling Under External Pressure of Ring Stiffened Circular 
Cylinders. In Transactions of The Royal Institution of Naval Architects VM1-15 vl07 139-156. 
United Kingdom: Royal Institution of Naval Architects. 


[9] Ross, C. T. F. 1965. The Collapse of Ring-Reinforced Cylinders Under Uniform External 
Pressure. In Transactions of The Royal Institution of Naval Architects VMM 5 vl07 375-394. 
United Kingdom: Royal Institution of Naval Architects. 

[10] Avallone, Eugene A. ed. and Theodore Baumeister III, ed. 1996. Marks ’ Standard 
Handbook for Mechanical Engineers. New York: McGraw-Hill. 

[11] Will, William. 2002. Private Conversations. Washington, D.C. 


57 


This page intentionally left blank 


58 


Appendix A: Codes for Test Cylinder l.d 
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AMERICAN BUREAU OF SHIPPING 


Rules for Building and Classing Underwater Vehicles, 


Systems, and Hyperbaric Facilities 

Definitions 


ksi := 1000|jsi 


E := 3000(ksi 

Modulus of Elasticity 

v := 0.3 

Poisson's Ratio 

Gy := 8000Q)si 

Yield strength 


Shell Parameters 

L s := 4.266n Distance between stiffeners 
L c = 22.488in Distance between bulkheads 

R = 8.007in Mean radius of shell 
t = 0.081 in Thickness of shell 

R„ := R + — Outer radius of shell 
° 2 

D 0 := 2 R 0 Outer Diameter of shell 

Ring Stiffeners 

t w := o.l 38n thickness of web of ring stiffener 

depth := 0.57in height of ring stiffener 

b := t w faying width of stiffener (from P&S for I beam stiffener) 

bf := o.Oin breadth of ring stiffener 

b 2 := bf - t w breadth of ring stiffener minus the web thickness 
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d := O.Oin 


ring stiffener flange thickness 


1 ^ : = l s - t w Unsupported spacing between stiffeners 

L:= maj ( L s- L b) 


L = 4.266in 


c, := .5- 


t w *depth + b 2 'd 
v t w -depth + b 2 d 


first centroidal height of ring stiffener 


Cj = 0.285in 

c 2 := depth - c, second centroidal height of ring stiffener 


C 2 = 0.285in 

h := cj - d distance from centroid of ring stiffener to nearest edge of flange 
Rj. := R + ,5-t + c 2 radius to centroid of ring stiffener 
Rj. = 8.332in 

A s := (t w -depth + b 2 -d) cross-sectional area of ring stiffener 


A s = 0.079in 2 


^ : = f-Y^bf-c j 3 - b 2 h 3 + t w -c 2 3 j moment of inertia of ring stiffener about its centroidal axis 


I r = 2.13x 10 3 in 4 


Rf := R + ,5-t + depth Radius to tip of the stiffener 
Rf = 8.617in 

Combined Plate and Ring Stiffener 


L e := min| 


1.5-Vra 
v °- 75L s ) 


effective length of shell plate 


61 




area of effective plate 


L e = 1.208in 

Ap := L t 

A sp := A P + A s area of plate and ring stiffener 

H c := depth + t height of combined plate and ring stiffener 

H c = 0.054ft 

B] := L - t w plate length minus the web thickness 

111 

t w H c 2 + B, t 2 + b 2 d (2 H c - d) 

c lc := 5 - 7 ;— - 77 - neutral axis of combined plate and ring stiffener 

t w -H c + Bjt + dfc^ 

from outer fiber of plate (R 0 ) 

Cj c = O.IOlin 

Rc := R ~ -5-t + c j c radius to centroidal axis of combined ring stiffener and shell 
Rc = 8.067in 

h : = 3 { Lc 1 c 3 - B r( c lc - ‘) 3 + bf '( H c - c lc) 3 - ( bf - ‘wlfHc - c lc“ d ) 3 _ 
moment of inertia of combined plate and shell 

I e = 0.009in 4 


General Equations 



M = 5.297 

1 
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e = 6.809 



Q = 3.405 

cosh (2- Q) - cos (2- Q) 

N - 

sinh(2-Q) + sin(2-Q) 

N = 0.997 

„ „ (sinh(Q) cos(Q) + cosh(Q) sin(Q)) 

sinh(2Q) + sin(2-Q) 

G = -0.081 


H ._ sinh(2Q) - sin(2 Q) 
sinh(2Q) + sin(2Q) 

H = 0.998 

Inter-Stiffener Strength (6.19.1) 

1) Inter-stiffener strength equations 


This equates to axisymmetric bucklim 



Effective area of plate and stiffener (External stiffeners) 


A = 0.073in 2 


F := 



A + t w -t + 


2N-t-L 

0 


F = -0.027 



yield pressure at midbay and midplane of cylinder 
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P y = 787.9psi 


This corresponds to Lobar bucklinr 


von Mises buckling pressure for a cylinder 


Uf. 


— -0.4$) — 
2 R l 2R 


P = 604.9psi 


Maximum allowable working pressure for inter-stiffener strength 


P P 

P c := -2 if < 1 
2 P 

y 

f ? \ p 

P. iM< > s3 

U 2P m p . 


5 r m 

- P y if->3 

6 y P y 


— = 0.768 

p y 


P c = 302.4psi 


P aits * *V 8 


P aits = 242 P si 
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2) longitudinal stress 

Limiting pressure corresponding to the longitudinal stress at stiffeners reaching yield. No direct 
correlation to major failure modes 


A I - - 

l 2 


A + t w -t + 


2-N-t-L 


2 'V 12 V 

P, :=-—• 1 + - y-H 

1 R 2 

1 - v 2 


P| = 732.7psi 


p als := p l'- 67 Maximum allowable working pressure for longitudinal stress 


P als = 4 90.88psi 


b Overall Buckling Strength (6.19.5) 


A 2 (n) := n - 1 


Aj(n) := 


-|(n 2 + xf 


Et EI e A 2( n ) 

Pnl( n > :=— A 1 ( n) +--- 

R R 3 L 
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Given 

n > 2 

n 2 := Minimiz^p n j, n) 
n 2 = 3.109 

n must be an integer 

n 2int :=round ( n 2>°) 


n 2int “ 3 

n 2flrst = n 2int - '®® d n 2int = ' 

n 2int - * otherwise 

n 2first = 2 

n 2first 
n 2int 
n 2int + 1 
n 2int + 2 

n 2prime, = 3 





5.127x 10 3 
1.522x 10 3 
1.974x 10 3 
3.022x 10 3 


psi 


p n := min(p n2 ) 
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p n = 1522.3psi 


P aol := Pn - 5 
P aol = 761.131psi 

Summary 

Axisymmetric Buckling 
P y = 787.9psi 
Lobar Buckling 
P m = 604.9psi 
General Instability 
p n = 1522.3psi 

n 2int = 3 



Germanischer Lloyd 

1988 Edition 

Chapter 2 - Submersibles 

MPa := lO^Pa ksi := lO^psi 

General Definitions 

E:= 3000(ksi 


L 1 := 4.26<5n 

stiffener spacing 

L 3 := 22.488n 

length of cylinder between bulkheads or lines of support 

s := .081in 

thickness of shell 

R := 8.00^n 

radius to centerline of shell 

v := 0.3 

Poisson's Ratio 

k := 80-ksi 

minimum yield stress of material 

Stiffener Dimensions 

:= O.Oin 

flange thickness 

d^-:= O.Oin 

width of flange from web to edge of flange 

b 2 := 2-d f 


d w := 0.57in 

Height of web 

e w := 0.138n 

web thickness 

b ;= e w 

width of stiffener ring in contact with shell 


A] := ef-(e w + 2 d f ) + e w d w cross-sectional area of stiffener ring 
Aj = 0.079in 2 
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H stiflf := d w + e f 


H stiff= 057in 


2 2 
2 e w H stiff + b 2' e f 

2 e w H stiff + b 2 e f 


C| = 0.285in 


c 2 : “ H stiff _C 1 


C 2 = 0.285in 


h 2 C 1 “ e f 


h 2 = 0.285in 


■)[( e w + 2-d f )-c, 3 - (2-d f )-h 2 3 + e w -c 2 3 


Centroidal Moment of inertia of ring stiffener 


Ij =2.13x 10 3 in 4 


e := c 2 + .5 s distance from stiffener centroid to center of shell 


e = 0.325in 


Effective Stiffener and Shell 


^etest ' b * -\Z2-R-s 


L e := L etest if L etest - L ] Effective length of shell (eqn 45a,b) 
Lj otherwise 


L e = 1.277in 


I Ii H - 

e A, 1 12 


Moment of Inertia of combined plate and shell 


I £ = 6.92x 10 3 in 4 
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H stiffc := H stiff + s Total height of stiffener and plate 


H stiffc = 0 - 65,in 


dj := s 


bs l • e w 


2 2 

1 e w‘^stiffc + b 2' e f + I ' s (^‘^stiffc “ s ) 

2 e w' H stiffc + b 2 e f +bs r s 


c j c = 0.47in 


c 2c * ^stiffc c lc 


c 2c = 0.181 in 


R 0 := R - .5 s + c 2c radius of stiffener ring centroid including effect of L e 


Rq = 8.148in 


A := A j- Modified area of stiffener ring 

R o 


A = 0.076in 


e l :=c 2c + * 5s 


ej = 0.222in 
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2.4 Asymmetric Buckling (Lobar but not named this) 



n := 2 initial guess 

s 2 (n 2 - 1 + Xj 2 ) 

12 R 2 -(l -v 2 ) 

' ' 2 2 
n — 1 + .5-Aj 

E-s-P nl (n) 

Pnl( n ) := --- 

Given 

n > 2 

:= Minimiz^p n j ,n) 

n a = 10.612 n must be an integer 
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n 2 :=Minimiz<(p nla> n) 





2.6 Symmetric Buckling (axisymmetric) 

More Definitions 

For Elastic-Plastic Region (Eqns are in the program) 

E s := E Secant Modulus = Young's Modulus for elastic region 

Ej := E Tangent Modulus = Young's Modulus for elastic region 

v p := v Poissons Ratio, elastic-plastic; = Poisson's ratio in elastic 
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a = 0.063- 


C 5 := a Lj calculation factor for symmetric buckling 


C 5 = 6.809 


2-s 2 E 

R 2 -i3-(l - v 2 ) 


critical pressure, elastic, calculation factor 


Past = 3716.21 lpsi 

C 5 := a Lj calculation factor for symmetric buckling 
C 5 = 6.809 


C 6 :=-VT^G 


c 7 := —>/1 + G 
' 2 


C 8 C 5 C 6 


C 9 c 5 c 7 


cosh( 

c s) : 

2 / 
- cost 

c 9 ) 

2 

coshl 


|sinh 

( C 8 

) cos 

( C 9 

)sin| 



4 , c 7 

cosh (Cg)sin(Cgj sinh(Cg)cos (qj 

C 7 _ C 6 

cosh (Cg) • s i nh (c 8 ) cos (c 9 ) • sin ( 4 ) 
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cos ■ sin^Q^ cosh (Cg)- sinh (Cg) 


Fo := 


2 cosh ^Cg)-sinh(Cg) cos^Q^sin^Q^ 


1 - v 


cosh ^Cg) • sin^C^ sinh (Cg) ■ cos (C9) 


F 4 : = 


^ 2 cosh (Cg)- sinh(Cg) cos sin(C9) 


C 10 :_ ' 


v^| 
1 - - 
2J 


s-Lj 


+ — + 


sL l F 1 


'i-±' 

V hj 


C H :- 


0.91 


1 - v 


°o : =' 


-pR 


°x : - a o’ 2 + C 10’ C 1 l’ F 4 


^ : -°o ( ! ~ C 10' F 2 + v C 10' C 11 F 4) 


°xi : -°o' 2 _c io c ir F 4 


ct x =-3.286x 10 Pa 


° 4 »i : = CT o-( 


1 “ C 10* F 2 “ v C 10' C l r F 4j 


Ox = -6.65lx IQ 8 Pa 


°<i> 

K o : =- 

CT x 


.8k = 64 ksi 


Determination factor (if>.8k then must use E t , I 


: s , etc) 
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Oj-=J ° x 2 + ®* 2 -°x°* 


cjj= 83.538ksi 


0 := " K o °x 


a = 96.5ksi 


e, = 3.215x 10 


o := k- .8 + .2tanh 5— e, - 4 
V k 1 


a = 79456psi 


k I fa 

e := -■ 0.8+ 0.2-atanh 5— 

E V Ik 


e = 3.215x 10 


k IE 

E s := —• 0.8+ 0.2- tanh 5— e 
e _ V k 


E s = 24548.3 lksi 


E. := E- 1 - tanh 5— e - 4 

V Ik 


E, = 1019.968ksi 


1 f 1 ^ 

V :=--V — 

p 2 l 2 J E 


v p = 0.336 





Paa : = 





trial2:= + 0 

tria!2= 1 


iterate on m until trial is < trial2 



Summary 

Axisymmetric Buckling 
p m = 923.8psi 

Lobar Buckling 
p cr = 605.9psi 

General Instability 
p g = 1244psi 
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Society of Naval Architects and Marine Engineers 


Submersible Vehicle Systems Design (1990) 
ksi := lOOC^si 

Input Section 


Lf:= 4.266n 

Length between frames 

Ljj := 22.488n 

Length between bulkheads 

R := 8.00^n 

Radius of cylinder to centerline of shell 

t := 0.08lin 

shell thickness 

Rq ;= R + .5-t 

radius to outside of shell 

R o = 8 

.047in 

E := 30000ksi 

modulus of elasticity 

Gy := 80-ksi 

minimum yield stress of material 

H := 0.3 

Poisson's ratio 

Ring Stiffeners 

t w := 0.13Sn 

w 

thickness of web of ring stiffener 

depth 0.57in 

total height of ring stiffener 

bf := O.Oin 

breadth of ring stiffener 

d := O.Oin 

ring stiffener flange thickness 


b 2 := bf - t w breadth of ring stiffener minus the web thickness 
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first centroidal height of ring stiffener 


2 2 
t w -depth + b 2 d 

t w -depth + bjd 
Cj = 0.285in 

c 2 := depth - C) second centroidal height of ring stiffener 

c 2 = 0.285in 

h := c i ~ d distance from centroid of ring stiffener to nearest edge of flange 

h = 0.285in 

R,. := R + .5-t + c 2 radius to centroid of ring stiffener 

Rj. = 8.332in 

A r := (*w de P th + ^2' d ) cross-sectional area of ring stiffener 

A r = 0.079in 2 

lr:= ( 3} '( bf ' c » _ b 2' h3 + V c 2 3 ) moment of inertia of ring stiffener about its centroidal axis 

Ir = 2.13x I0~ 3 in 4 

Combined Plate and Ring Stiffener 

A p ; = Lf t area of effective plate 

A sp := A p + A r area °f plate and ring stiffener 

H c := depth + t height of combined plate and ring stiffener 

B l := L f " l w P* ate length minus the web thickness 
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neutral axis of combined plate and ring stiffener 


V H c 2 + B r t 2 + b 2 -d-(2-H c -d) 

ClC '~ L V H c + B ] t + d b 2 

from outer fiber of plate (RO) 

Cj c = O.lOlin 

R c := R - ,5 t + c lc radius to centroidal axis of combined ring stiffener and shell 
R^ = 8.067in 

>e }[ L f c lc 3 * B l (=lc - *)’ + bf K - Clef - ( bf - >w) ( H c - '1C - d f] 
moment of inertia of combined plate and shell 
l e = 0.009in 4 

Axisymmetric yielding 



e = 6.809 

cosh ( 9 ) - cos(0) 
sinh(0) + sin(©) 

N = 0.997 



p = 1.126 
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Py = 781.382psi 


Lobar Buckling 



= 604.9psi 


General Instability 

rc*R 
m :=- 

Lb 

n := 1 



Given 

n > 1 
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n 2 := Minimiz^p cr j,n) 

n 2 = 3.109 n must be an integer 

n 2int :=round ( n 2’ 0 ) 
n 2int = 2 

n 2first := n 2int' 100 if n 2int= 1 n 2flrst = 2 

n 2 j nt - 1 otherwise 

n 2first 
n 2int 
n 2int + 1 
n 2int + 2 

Pcrl^ n 2primej^) 

Pcrl^ n 2primej^ 

Pcrl^ n 2prime 2 ^ 

Pcrl^ n 2prime^ 

p C r : = min (Pcr2) 

p cr = 1522.3psi 

Summary 

Axisymmetric Buckling 
Py = 781.4psi 
Lobar Buckling 
Pk = 604.9psi 
General Instability 
p cr = 1522.3psi 
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MIT 13A Professional Summer Submarine Design Trends 

PROGRAM TO COMPUTE SUITABILITY OF SUBMARINE DESIGN PARAMETERS. 

Ref: "Hull Material Trade Off Study", D Fox, Jan 94 
ksi h lOOtysi 

Global Variable Inputs: 

e s 0.0-in Eccentricity 

Material: 


a v = 8000Q)si 

yield stress of material 

E = 3000(ksi 

young’s modululs 

v = 0.3 

poisson’s ratio 

Geometry: 


R £ 8.003n 

shell radius 

Dh 2 R 

shell diameter 

Lf = 4.266n 

frame spacing 

Ls = 22.488n 

bulkhead spacing 

t p s 0.08 lin 

shell thickness 

tf = O.Oin 

flange tickness 

Wf = O.Oin 

flange width 

t w = 0.138n 

web thickness 

h w £ 0.57in 

web height 
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Compute areas: 


r, := r + — Frame Radius 

2 . 

Rj = 8.047in 

A f := tf• w f Frame flange, web area 

:= t\v'h\y 

A := Af + A w Frame Area = Flange + Web 
A = 0.079in 2 

PARTI SHELL YIELDING 

Von Sanden and Gunther (1952) 

PNA Section 8.4 

t\v’tp 

B :=-— 

A + tw t P Area ratio 


B = 0.124 


0 := Lf 


3.(l-v 2 ) 

(Rtp) 2 


Slenderness Parameter 


9 =6.809 


cosh(0) - cos(0) 

sinh(0) + sm(o) Deflection coefficient 


P:=- 


2-N 

r i n 

0.25 r 


A + t w tp 

i 

CO 

1 

< 

to 

1_ 

. 


(3 = 1.126 


Frame flexibility parameter 
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Frame deflection parameter 


r = 0.341 


•J e ( 9 ) , ( 0 ] (0 

smh — cos — + cosh — -sin — 

U; U; UJ U 


sinh(0) + sin(0) 


Bending effect (mem) 


H m =0.081 


H E :=-2 


, ( 0 I ( 0 ) . (e 

3 1 wj V2 J I 2 J 12 


sinh(0) + sin(0) 


Bending effect (bend) 


H e = 0.085 


P crl := f y o y 

P n <r- 1 psi 


sinh(o) - sin(o) 
sinh(o) + sin(o) 


<— 1000 psi 
P2 4— 2000 psi 

delta <- 5 psi 
limit <r~ 1 psi 
conv<- 1 psi 
1^0 

while j < 20 
I for i € 0 .. 2 


-- 1 1 + r,H M + v-H e i] 

-PR 

"7—[1 + v h£ i] 

1 rt 


~P. R 

& xxso *- 1 0.5 + F H £i 


Bending efffect near frame 


-P.R 

a xxsi*->0.5 - P H £i 

l P 

- p i R [ r ( 3 v 5 ' 

-1 -r- 1 + v - k 

tp l L l.-v 2 J . 

- p , R [ [ ( 3 \ 05 1 

-1 - r 1 - v - k 

tp . ll-v 2 j 
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Perl _ 


801.682 

14 


PeritPcrl 0 ' lP s * 

Peril = 801 7psi 

PART 2 LOBAR BUCKLING 

Windenberg Approx of Von Mises (1933) 
Assumes n lobes = Pi*D/L 
Collapse pressure: 



p cLB = 604.9psi 

PART 3 GENERAL INSTABILITY 

Corrected Bryant Formula (1954) for better model test correlation 
Pressure loading is: 

P:= pgDfSFgj 
P = 6697.99psi 

Compute effective frame spacing: 



y = 1.802 

Compute clear length: 
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LcLf t w 
nj := 0.5 yj 1 - y 

n\ = 0.4481 
Web thickness: 
n 2 := 0.5-/7T”y 

n 2 = 0.837 

2 2 

4 cosh(n]*9) -cos (n 2 -e) | 

0 cosh(n]-0)sinh(n]-0) cos(n 2 - 0 )-sin(n 2 - 0 ) 

"I n 2 

Fj = 0.23 
must be less 
than 1.00 

FI is almost a linear decreasing function for pressures from 1 to 2000 psi with an average value 
of 0.27. This will be used in the following analysis as the pressure is the unknown and therefore 
the above equations cannot be directly used. 

F! := 0.27 

Lgff : = l c F] + t w Effective shell plate length: 

L e ff = 1.253in 

Theoretical critical lobe number values are: 


i := 0.. 2 



Circumferential Lobes 


A e ff := L e fftp Effective plate area 
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A e ff = O.lOlin^ 


m:=7x— Longitudinal Lobes 


m = 1.119 


h w + tf] h w + t p A 

2 ' Af - Acff 

j / j / 

A e ff + A w + A f 


Frame-plate neutral axis (ref web centre+ toward 


flange): 


y na = - 0 . 183 m 


Uses Parallel Axis Theorm: Icor = I + Ad2 


Moments of inertia for plate,flange, web: 


1 _ I A l P + 

J pcor •- + A e ff- “ + y na 

V 2 / 


Kvcor ■— I\v + A w -(y na ) 


I _ T . A + 

ncor •“ k + Af* " y na 

V 2 


Total: 


90 



^eff Ipcor + hvcor + Ifcor 


. 4 

I e ff = 0.007in 

The critical Elastic General Instability pressure is: 



Min Pressure: 

f 3 

5.035x 10 

p cGl= 1.278x 10 3 

v 1.516x 10 3 

PcGI : = m'KPcGl) 

p cGI = 1278-lpsi 

Summary 

Axisymmetric Buckling 
Peril = 801 7psi 
Lobar Buckling 
p cLB = 604.9psi 
General Instability 
p cGI = 1278.lpsi 
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Analytic Solution 


Definitions 


ksi := 6.89475710 6 Pa 

E:= 3000Cksi 

Modulus of Elasticity 

R := 8.001n 

Radius of cylinder to centerline of shell 

D:= 2 R 

Diameter of cylinder 

Oy := 80ksi 

Yield strength 

L ;= 4.266n 

Length of supported cylinder 

Lh := 22.488n 

Distance between bulkheads 

M := -3 

Poison's ratio for Fe/Steel 

t := 0.08 lin 

Shell thickness 

Ring Stiffeners 


t w := O.J3Sn thickness of web of ring stiffener 
H := o.57in height of ring stiffener 

b := t w faying width of stiffener (from P&S for I beam stiffener) 

bf := o.Oin breadth of ring stiffener 

b 2 := bf - t w breadth of ring stiffener minus the web thickness 

d := O.Oin ring stiffener flange thickness 

Lbnash : = L - bf distance from flange edges 
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yH 2 + b 2 d2 

+ b 2 *d 


first centroidal height of ring stiffener 


Cj = 0.285in 


c 2 H " C 1 


second centroidal height of ring stiffener 


C 2 - 0.285in 


h := Ci - d 


distance from centroid of ring stiffener to nearest edge of flange 


h = 0.285in 


Rj. := R + ,5-t + c 2 radius to centroid of ring stiffener 


Rj. = 8.332in 


A r := (t w H + b 2 d) cross-sectional area of ring stiffener 


A r = 0.079in 


'X 'X *x 

I r := - !■( bf-Cj - b 2 -h + t w -c 2 1 moment of inertia of ring stiffener about its centroidal axis 


Ij. = 2.13x 10 3 in 4 


A eff :=A r ~ 
I R,. 


effective area of stiffener eqn [24a] from P&S 


A e ff = 0.076in 


a :=- ratio of effective frame area to shell area eqn [62] P&S 

L-t 


a =0.219 


p := — ratio of faying width to frame spacing eqn [62] P&S 


P = 0.032 
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Combined Plate and Ring Stiffener 

A p := L-t area of effective plate 

A sp := A p + A r area of plate and ring stiffener 

H c := H + t height of combined plate and ring stiffener 

H c = 0.054ft 

Bj := L - t w plate length minus the web thickness 

t w H c 2 + B,t 2 + b 2 d(2H c -d) 

c ]c := 5- neutral axis of combined plate and ring stiffener 

t w * H c + Bjt + dt>2 

from outer fiber of plate (Ro) 

c j c = O.lOlin 

R c := R + .5-t + c lc radius to centroidal axis of combined ring stiffener and shell 
R c = 8.148in 

*e : = ~[ L c lc 3 - B l( c lc " O' + bf ( H c " c, c ) 3 - (bf - t w ) (H c - c, c - d) 3 ] 
moment of inertia of combined plate and shell 
I e = 0.01 iif* 

Buckling of unreinforced shells 

von Mises bucking pressure: 

n := 1 

guess at number of waves around the circumference 
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Axisymmetric Buckling 


p cab ■ 


y <- 0 

1 im it 4 — I 0 p si 
t est <— 0 p si 
conv <— 1 p si 

j <- 0 

wh i le j :< 2 0 

l 


n i *- yV> - y 


i 


1 ; *- — V> + r 

e <- ^3 (l - n 2 )~r 


cosh (r| ]0) 2 — c o s ^ r| ^g) 2 


0 cosh^ri j-0^ sinh^r| j-0^ cos|ii 2 ®)' s > n (*1 2'®) 

I 11 2 

cosh (t] j-G^sin^ri 2®) sinh^ri j O^ cos^ri 2®) 

_^_2_ilj_ 

cosh^n ,- 0 )sinh(r| ,0) cos(r| 2 ' 6 )sin (r| 2 ©) 

Hi n 2 


< 1 3 

■n 2 

^ 1 

J 1 - U 2 '“"(l 

1 -0^-sin h ^rj | 0 ) 

c °s(r| 2 * 0 ) sin (rj 2 °) 


n 1 

^ 2 

cosh 

1 •0^-sin 2 ) 

sinh^r| |-0 ^cos(ti 2 *^) 

. , 1. 3 

n 2 

*1 1 

J 1 - u 2 cosh ( T ' 

j -0^-sin h ^r) j -0^ 

cos(ti 2 -0) sin (rj 2 e ) 


^ 1 

^ 2 


0-f> 


a + p + (l - p)-F j 

denom ) <— A -^F 2 + F 2 *P 4'0 — 2 -j.i)- 


fF^f) 

VN 1 “M / 


+ F 4 2 -(l-^-f|.i 2 


denom 2 


P c 2 


GM-—ra 

°>(r) 


Ji 


*4- denom 1 - denom 2 


break if |p c2 - test | < lim it 

*unrv&(iy 

.est <- p c2 
j <- j + I 

P c 2 


96 



( 885.4^ 

Pcab = l 2 J 

Paxiy - Pcab^'^P s ' 


Paxiy 


885.4psi 


Asymmetric Collapse (Lobar buckling) 

Windenberg Approx of Von Mises (1933) 
Assumes n lobes = Pi*D/L 



PcLB = 6049 P si 


General Instability of shells and rings 



Given 

n > 1 

n 2 Minimiz^p c Qj,n) 

n 2 - 3.109 n must be an integer 
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n 2int :=round ( n 2’°) 


n 2int - 3 


n 2pi 


rime* 


n 2int - 1 
n 2int 
n 2int + 1 
n 2int + 2 


p cgi2 : ~ 


J 

PcGI^ n 2prime^ 

PcGI^ n 2primej^ 

PcGI^ n 2prime,^ 

PcGl^ n 2prime^ 


Pcgibryant m ' n (Pcgi2) 


Pcgibryant = 1522x 10 P si 


Summary 
Lobar Buckling 

PcLB = 6049 P si 

Axisymmetric Buckling 

Paxiy = 885 - 4 P si 
General Instability 
Pcgibryant = 1522.3psi 
n 2int = 3 
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Appendix B: Codes for Test Cylinder l.f 
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AMERICAN BUREAU OF SHIPPING 


Definitions 

ksi := IOOQdsS 

e := 3000Cksi Modulus of Elasticity 

v := 0.3 Poisson's Ratio 

c y := 9850Q)si Yield strength 

Shell Parameters 

L s := 2.665n Distance between stiffeners 

L c := 42.129n Distance between bulkheads 

R := I7.3285in Mean radius of shell 
t := o.263in Thickness of shell 

R •= R + — Outer radius of shell 
° 2 

R 0 = 17.46in 

D 0 := 2 R 0 Outer Diameter of shell 

Ring Stiffeners 

t w := o.i 98n thickness of web of ring stiffener 
depth := i.025n height of ring stiffener 

b := t w faying width of stiffener (from P&S for I beam stiffener) 

bf := o.763n breadth of ring stiffener 

b 2 := bf - t w breadth of ring stiffener minus the web thickness 
d := ,263in ring stiffener flange thickness 
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L s l w 


Unsupported spacing between stiffeners 


Lj, = 2.467in 


L:= ma>(Ls,Lb) 


L = 2.665in 


c,:=.5. 


( ,2 L ,, 2 ^ 

t w * depth + b 2 *d 

t w - depth + b 2 *d 


first centroidal height of ring stiffener 


C| = 0351 in 


C 2 := depth - Cj 


second centroidal height of ring stiffener 


C 2 = 0.674in 

h := cj - d distance from centroid of ring stiffener to nearest edge of flange 

Rj := R - ,5-t - c 2 radius to centroid of ring stiffener 
Rj = 16.523in 

A s := (t w - depth + b 2 -d) cross-sectional area of ring stiffener 


A $ = 0.352in 2 

Ij. := Qj(bf.c, 3 - b 2 h 3 + t w -c 2 3 j moment of inertia of ring stiffener about its centroidal axis 
1,.= 0.031 in 4 

Rf := R - ,5-t - depth Radius to tip of the stiffener 
Rf = 16.172in 
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Combined Plate and Ring Stiffener 


L„ := mi 


l.sViTi 


Effective length of shell 


L e = 1,999in 


Ap := L t 


area of effective plate 


Asp ' Ap + A s 


area of plate and ring stiffener 


H c := depth + t height of combined plate and ring stiffener 


H c = 0.107ft 


B 1 ^ hv 


plate length minus the web thickness 


hv H c + Bj-t + b2’d-^2-H c - d^ 
t w *H c + Bjt + dt>2 


neutral axis of combined plate and ring stiffener 


from outer fiber of plate (R0) 


Cj c = 0.4in 


R c := R + .5 t - C| c radius to centroidal axis of combined ring stiffener and shell 


R c = 17.06in 


'« - j[ L -®|-( c lc- >) 5 + “ (He - «lc) 3 - (W “ 'w) ( H c - c,c - <l) 3 ] 


moment of inertia of combined plate and shell 


I e = 0.187in 4 


General Equations 


M = 1.248 
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1 



e = 1.605 


Q ; - 


0 ^ 

2 


Q = 0.802 

cosh(2*Q) - cos(2-Q) 
sinh(2*Q) + sin(2Q) 

N = 0.774 


^ (sinh(Q) cos(Q) + cosh(Q)-sin(Q)) 
sinh(2*Q) + sin(2Q) 

G= 0.934 

sinh(2Q) - sin(2 Q) 
sinh(2Q) + sin(2Q) 


H = 0.41 

Inter-Stiffener Strength (6.19.1) 

1) Inter-stiffener strength equations 
This equates to axisymmetric buckling 


A := A c 


f R ^ 




Effective area of plate and stiffener 


A = 0.369in 


2 



A- 


1 -- 
2 


•G 


A + t w -t + 


2N-1-L 

6 


F = 0.267 
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yield pressure at midbay and midplane of cylinder 


p := —-5. 

» 1-F 


Py = 2039.3psi 


von Mises buckling pressure for a cylinder 

This corresponds to Lobar buckling 


2.42 E- 


(.-vf- 


4 / \ 2 

■-0.4sf—1 

2 R l 2-R ) 


P m = 10369.9psi 


maximum allowable working pressure for inter-stiffener strength 


P P 

m r m 

- if -<1 

2 P 

y 

( p A p 

y m 

P v - 1-— if 1 <— <3 

y 2P P 

V ™) y 


5 r m 

—P if — > 3 

6 y p y 


= 5.085 


P c = 1699.4psi 


^aits • P c -8 


P ait , = 1359.5psi 


104 



2) longitudinal stress 

limiting pressure corresponding to the longitudinal stress at stiffeners reaching yield 
No direct correlation to major failure modes 


A + t w t + 


2-N-t-L 


2 a a 

P, :=-—• 1 + 

1 R 


Pi = 2098psi 


P , := P.-.67 maximum allowable working pressure for longitudinal stress 


P a]s = 1405.69psi 


Overall Buckling Strength (6.19.5) 


A2(n) := n - 1 


Aj(n) := 


A 2 (n) + — I- 


( 2 5i 2 ) 2 

•\n + A, / 

) 


Et EI eA 2 (n) 

p n i( n ) : = v' A i (n) + —;— 

R R 3 L 
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n > 2 

n 2 := Minimiz^p n j ,n) 

n 2 = 2.974 n must be an integer 

n 2int :=round ( n 2>°) 
n 2int = 3 


n 2firsf n 2int'^ d n 2int~ * 
n 2int “ ^ otherwise 



p n := min(p n2 ) 


p n = 4496. lpsi 
^aol := Pn 

P a0 l = 2.248x 10 3 psi 
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Summary 

Axisymmetric Buckling 
Py = 2039.3psi 
Lobar Buckling 
P m = 10369.Sfsi 
General Instability 
p n = 4496. lpsi 

n 2int = 3 
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Germanischer Lloyd 


MPa := 10 6 Pa 
ksi := 1 O^psi 

General Definitions 


E:= 3000Cksi 

modulus of elasticity 

L] := 2.665n 

stiffener spacing 

L 3 := 42.129n 

length of cylinder between bulkheads or lines of support 

s := 0.263in 

thickness of shell 

R := 17.3285n 

radius to centerline of shell 

v 0.3 

Poisson's Ratio 

k := 98.5ksi 

minimum yield stress of material 

Stiffener Dimensions 

e f :=0.263n 

flange thickness 

d f := 0.2825n 

width of flange from web to edge of flange 

bj * = 2 -dj- 


d w := 0.762n 

Height of web 

e w := 0.198(1 

web thickness 

b; =e w 

width of stiffener ring in contact with shell 


A l := e f( e w + 2 ' d f) + e w' d w cross-sectional area of stiffener ring 
A, = 0.352in 2 
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H stiff d w + e f 


H s tiff = 1 025in 


2 2 
1 e w H stiff + b 2 e f 

1 2 e w‘ H stiff + b 2 e f 

C| = 0.351 in 

c 2 := H stiff - c 1 

C 2 = 0.674in 


h 2 := c i - e f 


h 2 = 0.088in 


^ 3 


( e w + 2 d f)' c l 3 " ( 2 d f) h 2 3 + e w' c 2 3 _ 


Centroidal Moment of inertia of ring stiffener 


I] = 0.031 in 4 

e := c 2 + .5 s distance from stiffener centroid to center of shell 


e = 0.805in 

Effective Stiffener and Shell 


L etest b + V2-R-S 




Retest ^ Retest - ^1 
L| otherwise 


Effective length of shell (eqn 45a,b) 


L e = 2.665in 


A|-e 


1 + 


+ 1 1 + 


V s 


12 


L e -s 


Moment of Inertia of combined plate and shell 


I e = 0.187in 4 


109 



H stiffc := H stiff + s Total height of stiffener and plate 
^stiffc = 1 -288in 

d l :=s 

bs l :=L e- e w 

2 2 

1 e vv ^stiffc + ^2‘ e f + r s '(2'H s tiff c ~ s ) 

c ] c ■” ~~ -- 

" e w'^stiffc + ^2' e f + V s 

C| c = 0.888in 

c 2c :=H stiffc ~ c lc 
c 2c = 0.4 in 

R 0 := R + .5 s - c 2c radius of stiffener ring centroid including effect of L e 
R 0 = 17.06in 

R 2 

A := A | —- Modified area of stiffener ring 

R 0 

A = 0.363in 2 
e l :=c 2c + - 5 ' s 

2.4 Asymmetric Buckling (Lobar but not named this) 

TlR 

x \ '■= ~— eqn 19 

L l 

n := 2 initial guess 
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Pnl( n ) 


f 2 \ 

n 

-+ 1 

j 


s 2 (n 2 - 1 + X, 2 ) 2 
12 R 2 (] - v 2 ) 


2 

n 


- 1 + ,5-Xj 


2 


Pnl( n > : = 


E' s 'P n ](n) 


R 


Given 

n >2 

ng := Minimiz^p n | ,n) 


n = 10.936 n must be an integer 

o 


n gint : = round (n g> 0) 


n gint ” 11 


Ngfirst 


n gint ^ n gint * 
n gint - * otherwise 


n gfirst 


= 10 


gprime * 


f n gfirst " 


f 10 l 

n gint 


ii 

n gint + 1 

n gprime 

12 

n • + 2 


v 13. 

gmt T *■) 




Pn2 := 


^Pnl^ n gprime 0 j X 
Pnl^gprime^ 
Pnl ("gprime^ 


^ p nl^ n gprime 3 


Per := min(p n2 ) 
p cr = 9269. lpsi 


eqn 18 


eqn 17 


111 




3.3.3 Genera] Instability 

L := L 3 Spacing between effective stiffeners 

71 R 

V= T 

x 2 = 1 -292 

x 2 4 

P n2 (") : =-- 

(n 2 - I + 0.SA. 2 2 ^n 2 + X 2 ) 

Es.(3 n2 (n) 

(n 2 - l) E 1 
p,(n)- 

Pnla( n ) := P 0 ( n ) + P]( n ) 
n := 2 

Given 


n > 2 

n 2 :=Minimiz<(p n]a> n) 

n 2 = 2.955 n must be an integer 


n 2int := round(n 2 ,0) 
n 2int = 3 
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psi 




r 11569^ 
4651 
6612 
^10240 j 


2.6 Symmetric Buckling (axisymmetric) 

More Definitions 

For Elastic-Plastic Region (Eqns are in the program) 

E s = E Secant Modulus = Young's Modulus for elastic region 

Ej := E Tangent Modulus = Young’s Modulus for elastic region 

v D := v Poisson’s Ratio, elastic-plastic; = Poisson's ratio in elastic 

shape factor 

calculation factor for symmetric buckling 

critical pressure, elastic, calculation factor 

Past = 8364.87(psi 


1 2 \ 

HI - v j 


s 2 R 2 


a - 0.024- 


mm 


C5 := a-Lj 


C 5 = 1.605 


Past • 


2s 2 E 


R 2 -J3-(l -v 2 ) 
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G 


P 

Past 


c 7 -.= - 2 4TTg 

Cg := C 5 -C 6 
C 9 :=C 5' C 7 





cos(c 9 )sin(c 9 ) cosh(Cg)sinh(Cg) 


C 7 

C 6 

cosh (Cg)sinh(Cg) 

cos(c 9 )sin(c 9 ) 

C 6 

+ c 7 

cosh (Cg)'sin(c<^) 

sinh(Cg)*cos(c<^ 

C 7 

C 6 

cosh (Cg)sinh(Cg) 

cos^C^sin^Q^ 

C 6 

+ c, 




«x : -«o(^C ] 0 .C„.F 4 ) 

- C 10' F 2 + v C 10 C 11' F 4) 

0 x] := <v(j " C 1 O' C 1 l’ F 4^j 

a x = -8.275 x 10 8 Pa 

0 <|>1 := 0 o( ] _C 10 F 2 _ v C ]0' C 11 F 4) 

0 ^ = -1.013x 10 9 Pa 

if the calculated stress is > .8 of the yield stress, then must use E t , E s , etc 

.8 k = 78.8ksi 


ct 4> 

K 0 :=- 

0 x 


2 2 

°i := V CT x +<J 4> “ 0 x c 


a: = 135.505ksi 


a:=-K 0 a 


X 


a = 147ksi 


G 



Si = 4.898x 10 


g := k 


.8 + .2-tanh 



a = 98461. Ipsi 
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E, = 59.275ksi 



H 3»[('-2 v p)-(2- v pK] 

H 2 : =[(2-Vp)-('-2Vp).Ko] 

H 1»[' + H 4 [h 2 2 -3.(|-v p 2 )]_ 
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%% H 4 



Pm = 293L4 P si 


p = 293 lpsi 

iterate on p until p m is equal to p 
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Summary 

Axisymmetric Buckling 

p m = 293 , . 4 p si 

Lobar Buckling 
p cr = 9269. lpsi 
General Instability 
Pg = 4651.4psi 
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Society of Naval Architects and Marine Engineers 

Submersible Vehicle Systems Design 

ksi := 100ty>si 

Input Section 


L f :=2.665n 

Length between frames 

Lj, := 42.129n 

Length between bulkheads 

R := 17.3285n 

Radius of cylinder to centerline of shell 

t := 0.263n 

shell thickness 

Rq R + .5 t 

radius to outside of shell 

Rq= 17.46in 

E := 30000 ksi 

modulus of elasticity 

Gy := 98.5ksi 

minimum yield stress of material 

fi := 0.3 

Poisson's ratio 

Ring Stiffeners 

t... := 0.19Sn 

w 

thickness of web of ring stiffener 

depth := 1.025n 

total height of ring stiffener 

bf := 0.763n 

breadth of ring stiffener 

d := 0.263n 

ring stiffener flange thickness 

b 2 := b f - t w 

breadth of ring stiffener minus the web thickness 
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first centroidal height of ring stiffener 


f 22 

t w -depth + b 2 d 

c ] := -5- - 

^ t w -depth + b 2 d 

Cj = 0.351 in 

c 2 := depth - cj second centroidal height of ring stiffener 

c 2 = 0.674in 

h := C| - d distance from centroid of ring stiffener to nearest edge of flange 

h = 0.088in 

R,. := R - .5 t - c 2 radius to centroid of ring stiffener 
R,. = 16.523in 

A r := (t w depth + b 2 d) cross-sectional area of ring stiffener 
A r = 0.352in 2 

V = Q)( bf c, 3 - b 2 'h 3 + V c 2 3 ) moment of inertia of ring stiffener about its centroidal axis 
I,. = 0.03 Jin 4 

Combined Plate and Ring Stiffener 

A p := Lf-t area of effective plate 

A sp := A p + A r area P^ a t e an( l r * n g stiffener 

H c := depth + t height of combined plate and ring stiffener 

B, := L f - t w plate length minus the web thickness 
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c icM 


V H c 2 + B r t2 + b 2 d ( 2H c- d ) 
V H c + B r t + db 2 


neutral axis of combined plate and ring stiffener 


from outer fiber of plate (Ro) 


c lc = 0.4in 


R *= R + .5h - Ci 


radius to centroidal axis of combined ring stiffener and shell 


R c = 16.972in 


*e := j_ L T c lc - B 1 ( c lc - l ) +bf ( H c- c lc) -( bf - t wK H c- c lc- d )' 


moment of inertia of combined plate and shell 


I e = 0.187in 


a. Axisymmetric yielding 


A r + V 1 


B = 0.129 




2-R R 


0 = 1.605 


cosh(0) - cos(g) 
sinh(o) + sin(o) 


N = 0.774 
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H = -0.833 



b. Lobar Buckling 



p b = 10369.94£si 


c. General Instability 


m:= 


TtvR 


n := 1 



Given 

n > 1 

1*2 := Minimiz^p cr |, nj 

n 2 = 2.974 n must be an integer 

n 2int :=round ( n 2-°) 
n 2int = 3 

n 2first := n 2inf 100 if n 2int = 1 
n 2int “ ^ otherwise 

n 2first “ 2 





11510.622' 
4496.065 

psi 

6320.305 
9773.51 ) 


Per : = mitl (Pcr2) 
p cr = 4496. lpsi 
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Summary 

Axisymmetric Buckling 
Py = 1927.7psi 
Lobar Buckling 
Pjj = 10369.^>si 
General Instability 
p cr = 4496.Ipsi 
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MIT 13A Professional Summer Submarine Design Trends 

PROGRAM TO COMPUTE SUITABILITY OF SUBMARINE DESIGN PARAMETERS. 

Ref: "Hull Material Trade Off Study", D Fox, Jan 94 
Define input parameters: 
ksi = lOOQ^si 

Global Variable Inputs: 


e = 0.0-in 

Eccentricity 

Material: 

Gy = 9850Q)si 

Yield Strength 

E = 3000(ksi 

Modulus of Elasticity 

v = 0.3 

Poisson’s Ratio 

Geometry: 

D = 34.657in 

shell diameter 

D 

R = — 

2 

shell radius 

tf = .263in 

flange tickness 

Lf = 2.665n 

frame spacing 

wf = .763in 

flange width 

4 = 42.129in 

bulkhead spacing 

t w = . 198in 

web thickness 

tp = 0.263n 

shell thickness 

h w = .762in 

web height 
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Ip 

R| := R —— Frame Radius 

2 

Rl = 17.197in 

Af := tf- w f Frame flange, web area: 

A\v := bv'^w 

A := A f + A w Frame Area = Flange + Web 

A = 0.352in 2 


PARTI SHELL YIELDING 


Von Sanden and Gunther (1952) 
PNA Section 8.4 


B := ■ 


bv'tp 


A + t w t p 


Area ratio 



cosh(6) - cos (e) 
sin h(9) + sin(9) 


Slenderness parameter 


Deflection coefficient: 


P:=- 


2-N 


A + t w t p 


r i i 

0.25 r 

\ - 

UJ 

1 

< 

_ NJ 

1_ 

V 


Frame flexability parameter: 


P = 1.676 



l + p 


Frame deflection parameter: 


r = 0.269 
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sinh 


Bending effect(mem): 


Hm==-2- 



+ cosh 




V 


e 

•sint — 


sinh(9) + sin(0) 


H m = -0.934 


0.5 sinh 


H E := ~2* 


2 

1 - v j 


— -cos — 


cosh 


fi 

u > 


0 

sinl — 


sinh(o) + sin(o) 


H e = 0.368 


sinh(o) - sin(o) 
sinh(o) + sin(o) 


Bending efffect near frame: 


P crl •= 


f y ^ y 

P 0 1 psi 
Pj 1000 psi 
<r~ 5000 psi 


delta <— 5 psi 
limit 1 psi 
conv <— 1 psi 
j^O 

while j < 200 
I for i e 0.. 2 


a <t«J>so 

a xxso 

CT xxsi 

a 4*|ifo *“ 


G 4>ifi <~ 


-P.R 

i 

l P 

-P.R 

l 

"V 

-P. R 
1 

-P. R 
1 

"V 

-P.R 

1 

1 P 

-P.R 

"V 


[l + T!Hm + v-H 

[l + r-'H M - v-H E'J 


• i 0.5 + r-H£! 


(0.5 - r-H E i 


' 


f 3 \ 05 ]1 

i - r- 

1 + V 

1 3 1 K 





_ 

U - V ) JJ 

" 


' 3 5 T1 

i - r 

1 - V 

K 






,i-v; jj 


Bending effect (bend) 
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- p i R r / 3 -\° 5 i 

—r r r v (—) k J 

lp L L li-v 2 J J_ 


- r i^ / 3 \ U:> 

- 0 5 - r - - - I K 

^ L It - v^J 

- p i R r , 3 \ 0.5 

- 0 . 5 + r [ - K 

l * L U-v 2 J 


h ^° s y 0 


°3 CT 

CT 4 O S y 7 


a SYM<-V a l -°l C5 2 + V 2 


( 2 2 ^ 
a SYF l a 3 -03-04 + 04 , 


cj 5Y 4— max 

\\ a SY 
stress. 4— a 5Y 

test <— 1 f y - g SY 1 


break if jstress^ - fy| < limit 
Pq <— Pj if shel dav > 0 
?2 ^1 ^ shel dav < 0 


P 1 p o + 

j«-i+ 1 


* ?2 “ 


Out, 4- j ~ 1 
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1994.1 P 

16 j 


Peril : = Pcrl 0 IPS' 


Pc^ = 1994. lpsi 


PART 2 LOBAR BUCKLING 

Windenberg Approx of Von Mises (1933) 
Assumes n lobes = Pi*D/L 
Collapse pressure: 


p cLB := -— 

f 4 




2.42 E- 



.2.5 



p cLB = 10369.^)si 


PART 3 GENERAL INSTABILITY 

Corrected Bryant Formula (1954) for better model test correlation 
Pressure loading is: 

P := 600Cf)si 

Compute effective frame spacing: 



y = 0.717 

Compute clear length: 

Lc := Lf ~ t\v 
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n ] := 0.5 >71 - y 
n, = 0.266 

Web thickness: 
n 2 := 0.5-V 1 + y 


n2 = 0.655 
t w - 0.198in 


Fl 


£ 

9 


2 2 
cosh(nj-O) -cos (n 2 e) 

cosh (nj-O^sin h(n t e) COS (n 2 -0) sin(n2-0) 

n l n 2 


must be less than 1.00 


Fi = 0.958 


Fi is almost a linear decreasing function for pressures from 1 to 10000 psi with an average value 
of 0.91. This will be used in the following analysis as the pressure is the unknown and therefore 
the above equations cannot be directly used. 

F) := 0.96 


Leff : = 4 F i + »w Effective shell plate length: 
Leff = 2.566in 

Theoretical critical lobe number values are: 


i := 0.. 2 

f 2} 


n := 


3 

V4, 


Circumferential Lobes 


A eff := 4fftp Effective plate area: 


Aeff = 0.675in 


2 
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m := n — Longitudinal Lobes: 

Ls 


m= 1.292 


h\v + *f | [ h w + *p | 

-I- - A f- -I- A eff 

2 y 2 / 

A eff + A w + A f 


Frame-plate neutral axis (ref web centre+ toward flange): 


Yna= -0.237in 


Uses Parallel Axis Theorm: Icor = 1 + Ad2 


Moments of inertia for plate,flange,web: 


Ipcor *p + A e ff 


tp + h\v 


2 

hvcor := I\v + A w *(y na ) 


f tf + h w 

Ifcor If + Af* — Yrja 

V 2 


Total: 


leffIpcor + I\vcor + Ifcor 


Moment of Inertia corrected for neutral axis. 


leff = 0.185in 4 


The critical Elastic General Instability pressure is: 
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P cG i := min(P cG i) 


P C GI = 4460.2psi 
n gi :=3 

Summary 

Axisymmetric Buckling 

Pcrit = 1994. lpsi 
Lobar Buckling 
PcLB = 10369.^si 
General Instability 
P C GI = 4460.2psi 
n gi = 3 


132 


Analytic Solution 


Definitions 

ksi := 6.89475710 6 Pa 
E := 3000(ksi 

R := 17.3285n 

D := 2-R 
Gy 98.5ksi 
L := 2.665n 
Lh := 42.129n 
»:= 3 
t := 0.263in 
Ring Stiffeners 

t w :=0.19Sn 
H := 1.025n 
b: = t w 
bf := 0.763n 

b 2 : =bf-t w 

d := 0.263n 
Lbnash := L _ bf 


Modulus of Elasticity 

Radius of cylinder to centerline of shell 

Diameter of cylinder 

Yield strength 

Length between stiffener centers 
Distance between bulkheads 

Poison's ratio for Fe/Steel 
Shell thickness 

thickness of web of ring stiffener 
height of ring stiffener 

faying width of stiffener (from P&S for I beam stiffener) 
breadth of ring stiffener 

breadth of ring stiffener minus the web thickness 

ring stiffener flange thickness 
distance from flange edges 
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c, := 5 


t vv H 2 + b 2 d 2 

t w -H + b 2 d 


first centroidal height of ring stiffener 


Ci = 0.351 in 


c 2 := H - cj second centroidal height of ring stiffener 


c 2 = 0.674in 


h:=c, -d 


distance from centroid of ring stiffener to nearest edge of flange 


h = 0.088in 


R,. := R - ,5 t - c 2 radius to centroid of ring stiffener 


R r = 16.523in 


A r := (t w H + b 2 d) cross-sectional area of ring stiffener 


A r = 0.352in 2 


I 1 H kf 3 U U 3 3 

\ = “ | bf c l - b 2 b + V C 2 


moment of inertia of ring stiffener about its centroidal axis 


1,.= 0.031 in"^ 


A efT= A r - 


effective area of stiffener eqn [24a] from P&S 


A^ff = 0.369in 


ratio of effective frame area to shell area eqn [62] P&S 


a = 0.526 


P := - ratio of faying width to frame spacing eqn [62] P&S 


p = 0.074 
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Combined Plate and Ring Stiffener 

A p := L t area of effective plate 

A sp := A p + A r area of plate and ring stiffener 

H c := H + t height of combined plate and ring stiffener 

H c = 0.107ft 

Bj := L - t w plate length minus the web thickness 

t w H c 2 + B, t 2 + b 2 d (2 H c - d) 

Ci»:= .5- neutral axis of combined plate and ring stiffener 

1 V H c + B r t + db 2 

from outer fiber of plate (RO) 

c lc = 0.4in 

R c := R - .5 1 - c lc radius to centroidal axis of combined ring stiffener and shell 

R c = 16.797in 

!e : = }[ L c lc 3 " B l ( c lc - O' + bf ( H c - c lc ) 3 - ( bf - VK H c - c lc - d ) 3 ] 
moment of inertia of combined plate and shell 

I e = 0.187in 4 

a. Buckling of unreinforced shells 

von Mises bucking pressure: 

n := 1 guess at number of waves around the circumference 
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Purl*") 



1 




n + 


V l ; 



Given 


n > 1 

n j := Minimiz^p ur j ,n) 

nj = 10.906 n must be an integer 

"lint := round (n,,0) 

"lint = 11 



p ur := min(p url ) 


p ur = 9272.3psi 
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c. Asymmetric Collapse (Lobar buckling) 

Windenberg Approx of Von Mises (1933) 
Assumes n lobes = Pi*D/L 



PcLB “ 10369.9psi 


d. General Instability of shells and rings 



n := 1 

Given 

n > 1 

n 2 := Minimiz^p c Q|,n) 

n 2 = 2.974 n must be an integer 
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n 2 int : = round (n 2 ,0) 


1 2int : 


ri 2prirne’ 


/n 2int _ 1 ^ 


n 2int 
n 2int + 1 
n 2int + 2 


Pcgi2 : 


PcGI^ n 2prim^ 
PcGI^ n 2primej^ 
PcGl^ n 2prim^ 
PcG^ n 2prim^^) 

Pcgibryant m ' n (Pcgi2) 


Pcgibryant =4496 -'P si 


Summary 

Lobar Buckling 
PcLB = 10369.%>si 
Axi symmetric Buckling 

Paxiy = 2140 - 7 P si 
General Instability 
Pcgibryant = 44 96.1psi 
n 2int = 3 
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Appendix C: Codes for Test Cylinder 2.a 
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AMERICAN BUREAU OF SHIPPING 

Rules for Building and Classing Underwater Vehicles, Systems, and Hyperbaric Facilities 
Definitions 


ksi := lOOQjsi 


E:= 3000(ksi 
v := 0.3 
Oy := 6550(psi 

Shell Parameters 


Modulus of Elasticity 
Poisson's Ratio 
Yield strength 


L s := 1,366n 
L c := 8.636n 
R := 8.4179n 


t := 0.0858n 


R ° :=R+ i 


Distance between stiffeners 

Distance between bulkheads 

Mean radius of shell 
Thickness of shell 

Outer radius of shell 


R 0 = 8.461 in 

D o : = 2 Ro 
Ring Stiffeners 

t w := 0.044n 
depth := 0.5319n 
b:= t w 
bf := 0.399n 
b 2 :=bf-t w 


Outer Diameter of shell 

thickness of web of ring stiffener 
height of ring stiffener 

faying width of stiffener (from P&S for I beam stiffener) 
breadth of ring stiffener 

breadth of ring stiffener minus the web thickness 
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d := 0.078n 


ring stiffener flange thickness 


:= L g - t w Unsupported spacing between stiffeners 

L b = ].322in 
L:= ma^Lg.Lj 

L = 1.366in 


( 2 2 
t w *depth + 

v V'depth + b 2 d 


first centroidal height of ring stiffener 


Cj = 0.143in 


c 2 := depth - c j 


second centroidal height of ring stiffener 


c 2 = 0.389in 

h := cj - d distance from centroid of ring stiffener to nearest edge of flange 

Rg := R + ,5-t + c 2 radius to centroid of ring stiffener 
Rg = 8.85in 

A s := (t w -depth + b 2 -d) cross-sectional area of ring stiffener 
A s = 0.051in 2 

moment of inertia of ring stiffener about its centroidal axis 

1,.= 1.219x 10“ 3 in 4 

Rf := R + ,5-t + depth Radius to tip of the stiffener 

Rf = 8.993in 


W3 


3 , .3 ^ 3 

bf'C] - b 2 -h + V c 2 , 
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Combined Plate and Ring Stiffener 



L e = 1.025in 

A p := L t area of effective plate 

A sp : = A p + A s area of plate and ring stiffener 

H c := depth + t height of combined plate and ring stiffener 

H c = 0.051ft 

B i := L - t w P late length minus the web thickness 

, V H c 2 + B r t2 + b 2 d (2 H c - d) 

c ic : = 5 t H +B, t + d b 2 - neutral axis of combined plate and ring stiffener 

from outer fiber of plate (R 0 ) 

Cj c = 0.174in 

: = R - -5 t + c ]c radius to centroidal axis of combined ring stiffener and shell 
Rc = 8.549in 

:= }[ L c lc 3 " B l ( c lc “ l ) 3 + bf ( H c " c lc) 3 " ( bf - twK H c - c lc " d ) 3 . 
moment of inertia of combined plate and shell 

l e = o.oosin 4 

General Equations 

M :=— 

VrI 

M = 1.607 
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e = 2.066 



Q= 1.033 

cosh(2*Q) - cos(2-Q) 
sinh(2Q) + sin(2*Q) 

N = 0.942 


G , = 2 (sinh(Q) cos(Q) + cosh(Q) sin(Q)) 
sinh(2*Q) + sin(2-Q) 

G= 0.835 

n ._ sinh(2Q) - sin(2-Q) 
sinh(2Q) + sin(2Q) 


H = 0.631 

a Inter-Stiffener Strength (6.19.1) 

1) Inter-stiffener strength equations 
This equates to axisymmetric buckling 


A:=A S - 


( RV 


External stiffeners Effective area of plate and stiffener 


A = 0.046in 


F:= 


Ad 1-l*G 

2 


A + t w *t + 


2m 

e 


F = 0.209 
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yield pressure at midbay and midplane of cylinder 


P = 

y i _f 


Py = 844.1 psi 


This corresponds to Lobar buckling 

von Mises buckling pressure for a cylinder 


(,-v 2 ) 4 . -L-0.4-L 
2 R l 2-R 


P = 2947.7psi 


maximum allowable working pressure for inter-stiffener strength 


P m P m 

p := J5 if JH<, 

2 P 

y 

f p 'N p 

p t . 1--2- if I < — £ 3 
y l 2p n.J p y 

c P 

j ^ m 

— *P V if-> 3 

6 y P 

y 


= 3.492 


P c = 703.4psi 


P aits * 


P aits = 562.7psi 
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2) longitudinal stress 

limiting pressure corresponding to the longitudinal stress at stiffeners reaching yield 
No direct correlation to major failure modes 






Pj = 848.5psi 

maximum allowable working pressure for longitudinal stress 


P a | s = 568.47psi 


b Overall Buckling Strength (6.19.5) 
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Given 



Summary 

Axisymmetric Buckling 
P y = 844. lpsi 
Lobar Buckling 
P m = 2947.7psi 

General Instability 
p n = 6498.2psi 

n 2int = 4 
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Germanischer Lloyd 


1988 Edition 

Chapter 2 - Submersibles 

MPa := 10 6 Pa 
ksi := lO^psi 

General Definitions 


E:= 3000Qcsi 
L, := 1,366n 

L 3 := 8.636n 

s := 0.0858n 

R := 8.4179n 

v := 0.3 

k := 65.5ksi 

Stiffener Dimensions 


Modulus of Elasticity 
stiffener spacing 

length of cylinder between bulkheads or lines of support 

thickness of shell 
radius to centerline of shell 
Poisson’s Ratio 

minimum yield stress of material 


e^-:= 0.078n 
d f := ,1775n 
b 2 := 2 d f 
d w := 0.4539n 


e w := 0.044n 


flange thickness 

width of flange from web to edge of flange 

Height of web 
web thickness 


b := e 


W 


width of stiffener ring in contact with shell 


150 




cross-sectional area of stiffener ring 


A l := e f'( e w + 2 d f) + e w' d w 
2 

A j = 0.051in 

H stiff := d w + e f 

H s ti ff = 0.532in 

_2 e w H stiff 2 + b 2 e f 2 
1 2 e w‘ H stiff + b 2 e f 

Cj = 0.143in 

c 2 :=H stiff ~ C 1 

C 2 = 0.389in 

h 2 : = c , - e f 

h 2 = 0.065in 

ij : = j-[(e w + 2d f ).cj 3 - (2-d f )-h 2 3 + e w -c 2 3 J Centroidal Moment of inertia of ring stiffener 
Ij = 1.219X 10 _ 3 in 4 

e := c 2 + .5 s distance from stiffener centroid to center of shell 

e = 0.432in 

Effective Stiffener and Shell 

^etest := ^ + V2 R S 



31,645mm 


L e := L etest if L etest - L i Effective length of shell (eqn 45a,b) 

L| otherwise 

L e = 1.246in 
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+ 1 1 + 


Moment of Inertia of combined plate and shell 


I e = 7.732 x 10 3 in 4 


^stiffc ' ^stifif + s 


Total height of stiffener and plate 


H stiffc = 0 - 618in 


d,:=s 


bs 1 := L e ~ e w 


1 e vv H stiffc + b 2' e f + 


bs r s ( 2H stiffc ~ s ) 


e w ^stiffc + b 2’ e f + bs I s 


c j c = 0.435in 


c 2c * ^stiffc c lc 


c 2c ~ 0.183in 


Rq := R - 5 s + c 2c radius of stiffener ring centroid including effect of L € 


Rq = 8.558in 


A := A, - 


Modified area of stiffener ring 


A = 0.049in 2 


e l c 2c + - 5 s 


ej = 0.225in 



n > 2 

n g :=Minimiz^p nl ,n) 

n g = 13.832 n must be an integer 
n gint := round (n g ,o) 
n gint = ^ 

n gfirst ' = n gint‘^® n gint _ ^ 
n gint ” * otherwise 


n gfirst ^ 


n gprime *' 


n gfirst 
n gint 
n gint + 1 
n gint + 2 
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3.3.3 General Instability 

L := Lj Spacing between effective stiffeners 

, n-R 

X 2 = 3.062 

p n2 (n) :=-- 

^n 2 - 1 + 0.5X 2 2 )-(n 2 + 

Esp n2 (n) 

p o (n): =-^- 

(n 2 - l) E I e 

Pl (n):=- - - 

R o L] 

PnlaW := P 0 ( n ) + p l( n ) 
n := 2 


Given 


r >2 := Minimiz<^p n ] a ,n) 

n 2 = 4.01 n must be an integer 

n 2int :=round ( n 2>°) 
n 2int = 4 

n 2int" 1 
n 2int 
n 2int + 1 
n 2int + 2 

Pnla^ n 2prim^y 

Pnla^ n 2primej^ 

Pnla^ n 2prime,^ 

Pnla^ n 2prime^ 

Pg :=min (Pg2) 

Pg = 6184.9psi 

2.6 Symmetric Buckling (axisymmetric) 

More Definitions 
Default values for Elastic-Plastic Region (Eqns are below) 

E s := E Secant Modulus = Young's Modulus for elastic region 

E t := E Tangent Modulus = Young's Modulus for elastic region 

v := v Poisson’s Ratio, elastic-plastic; = Poisson's ratio in elastic 
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1 3 Vl - v 


shape factor 


a = 0.06- 


Cj-a-Lj 


calculation factor for symmetric buckling 


C 5 = 2.066 


2-s 2 E 

R 2 J3-(l - v 2 


critical pressure, elastic, calculation factor 


Past ~ 3772.56’^>si 

Must guess pressure and then change when compared to p m at bottom of calculation 


C 6 := T'V 1 ~G 


Cy — -yj 1 + G 


C 8 C 5 C 6 


C^:- C$Cj 


cosh( 

C 8> 

2 - COS ( 

c 9 ) 

2 

cosh j 

. C 8) 

|-sinh 

( C 8 

;) cos 


)sin| 

N 


H c 7 

cosh(Cg)sin(c^ sinh(Cg)cos^Cp) 

°7 _ C 6 

cosh (Cg)sinh(Cg) cos^C^jsin^C^j 
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cos (c<^sin(c<^ cosh (Cg)- sinh^Cg) 


F 3== 


^ 2 cosh ^Cg^sinh^Cg^ cos{c^j sin(Cgj 


cosh (Cg) * sin sinh (c^j ■ cos 


F 4 := 


^ 2 cosh^Cg)-sinh(Cg) cos^Q^sin^C^ 


C 10 : “' 


( A 

1 -- 

^ 2) 


s-Li 


-+ — + 

s-L| Lj 


v'-^r 


Cjl :- 


0.91 


CT o := 


1 -v 
p-R 


CT x : -°o'| 2 + C 10' C ir F 4 


a x =-4.675 x 10 Pa 


_C 10' F 2 + vC 10 C ir F 4) 


cr^ = -5.789x 10 Pa 


°xi : =<v _ c io‘ c i r F 4 


o xl = -2.292x 10 Pa 


CT <t>iM 1 " c iff F 2- v * c i(r c ir F 4) 


CTxj = -5.074x 10 Pa 
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,8k = 52.4ksi 

If the calculated stress is greater than ,8k then must find E t , E s , etc. 



e, = 2.799 x 10 



E, = 480.6ksi 



v p = 0.345 


158 




Cj = 0.555 


c r c 2~ v p c 3 2 


c 0 = 0.163 
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71 


trial = 0.353 

Iterate on m until trial <= trial2 


m:= 1 

- frn 
tnal2:= I— (m + 1) 

trial2 = 1 



Summary 

Axisymmetric Buckling 
p m = 1029.8psi 
Lobar Buckling 
p cr = 2814.7psi 
General Instability 
p g = 6184.9psi 
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Society of Naval Architects and Marine Engineers 


Submersible Vehicle 

Systems Design 

ksi := lOOQisi 


Input Section 


L f := 1.366n 

Length between frames 

8.636n 

Length between bulkheads 

R -8.417911 

Radius of cylinder to centerline of shell 

t := 0.0858n 

shell thickness 

Rq *— R + .54 

radius to outside of shell 


Rq = 8.461in 


E := 30000 ksi 
Gy := 65.5-ksi 
|i := 0.3 

Ring Stiffeners 

t w := 0.044n 

w 

depth := 0.5319n 
bf := 0.399n 

d 0.07Sn 
b 2 := bf - t w 


modulus of elasticity 
minimum yield stress of material 

Poisson's ratio 

thickness of web of ring stiffener 

total height of ring stiffener 

breadth of ring stiffener 

ring stiffener flange thickness 

breadth of ring stiffener minus the web thickness 
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first centroidal height of ring stiffener 


2 2 
t w - depth + b 2 ’d 

t w depth + b 2 'd 
Cj = 0.143in 

c 2 := depth - c j second centroidal height of ring stiffener 

C2 = 0.389in 

h := c, - d distance from centroid of ring stiffener to nearest edge of flange 

h = 0.065in 

R,. := R + .5-t + c 2 radius to centroid of ring stiffener 
Rj. = 8.85in 

A r := (t w depth + b 2 d) cross-sectional area of ring stiffener 
A r = 0.051 in 2 

j--(bf-cj 3 - b 2 h 3 + t w -c 2 3 j moment of inertia of ring stiffener about its centroidal axis 

1,.= 1.219X 10“ 3 in 4 

Combined Plate and Ring Stiffener 

A p := Lj-t area of effective plate 

A sp := A p + A r area °f P* ate anc ^ r i n 8 stiffener 

H c := depth + t height of combined plate and ring stiffener 

Bj := Lf - t w plate length minus the web thickness 
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neutral axis of combined plate and ring stiffener 
from outer fiber of plate (Ro) 

C] c = 0.174in 

: = r - , 5 -t + c lc radius to centroidal axis of combined ring stiffener and shell 
R c = 8.549in 

'e - - B l( c lc - <) 3 + bf ( H c - c,c) 3 - (bf - <w) fc “ c lc “ d ) 3 ] 

moment of inertia of combined plate and shell 
I e = 0.008in 4 


Axisymmetric yielding 



9 = 2.066 

cosh(0) - cos(g) 
sinh(0) + sin(0) 

N = 0.942 


V H c 2 + B r t 2 + b 2 -d-(2.H c -d) 

V H c + B r l + db 2 
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Py = 814.733psi 


Lobar Buckling 



P(j = 2.948x 10 3 psi 

General Instability 

tiR 
m:=- 

Lb 
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E-t 


m 


( 

2 ^ / v 

2 

, m ( 2 2) 

n - 

1 H-\n + m / 

V 

2 ) 


Pcrl( n > := 


Given 

n > 1 

n 2 := Minimiz^p cr j ,n) 

n 2 = 3.957 n must be an integer 

n 2int :=round ( n 2’°) 


(n 2 - l)-E-I 


n 2int = 4 


n 2 first ■ 


n 2int' 1°° if n 2int= 1 
n 2int “ * otherwise 


1 2prime' 


r n 2first ^ 


n 2int 
n 2int + 1 
V n 2int +2 y 


Pcr2 * 


Pcrl^ n 2prime^) 

Pcrl( n 2prime,) 


Pcrl( n 2prim^) 
Pcrl^ n 2prime^ 

Per :=min (Pcr2) 

p cr = 6498.2psi 
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Summary 

Axisymmetric Buckling 
Py = 814.7psi 
Lobar Buckling 


Pj, = 2947.7psi 


General Instability 
p cr = 6498.2psi 
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MIT 13A Professional Summer Submarine Design Trends 

PROGRAM TO COMPUTE SUITABILITY OF SUBMARINE DESIGN PARAMETERS. 

Ref: "Hull Material Trade Off Study", D Fox, Jan 94 
Define input parameters: 
ksi = 100C£>si 

Global Variable Inputs: 

e = o.o-in Eccentricity 

Material: 

cr y s 6550Qjsi Yield strength 

E s 3000(ksi Modulus of Elasticity 

v - °- 3 Poisson’s Ratio 

Geometry: 

R s 8.4179n shell radius 

D = 2 R shell diameter 

4 s i.36fin frame spacing 

4 s 8.63fin bulkhead spacing 

t p = o.0858n shell thickness 

tf = o.07Sn flange tickness 

Wf = o.399n flange width 

t w = 0.044n web thickness 

h w = o.4539n web height 
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Compute areas: 

t p 

Rf := R + — Frame Radius 

2 

Rf = 8.461 in 

Frame flange, web area: 

Af := tf- Wf 

A\v := t\v'h\v 

A := A f + a vv Frame Area = Flange + Web 

A = 0.051 in 2 


PART 1 SHELL YIELDING 



Area ratio 


0 := Lf- 



0 = 2.066 


cosh(Q) - cos (e) 
sinh(o) + sin(o) 


Slenderness parameter: 


Deflection coefficient: 


Z- N 1 13 

P := ~ 7 ~l - rr V R t p Frame flexibility parameter: 

A + V t p[ 3 .(i_v 2 JJ 

(3 = 1.947 
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V 

1 - - 
2 


-B 


Frame deflection parameter: 


r := 


1 + p 


r = 0.265 


sinh 


H M :=-2- 


H m = -0.835 


H E := -2- 




cos 






., e) . f o' 

+ cosh | — I-sin — 

2 Bending effect (mem): 


sinh(0) + sin(o) 


,0.5 sinh 


1 -v 


fe 

— |cos 

u 


— - cosh 

V2 


. fe 

\2j ^ 2 


sinh(9) + sin(9) 


Bending effect (bend): 


H e = 0.557 

sinh(o) - sin(o) 
sinh(o) + sin(o) 


Bending efffect near frame: 
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P crl : ~ 


f y <r~ a y 

Pq «- 1 psi 
Pj 1000 psi 


<— 2000 psi 

delta <r- 5 psi 
limit <r~ 1 psi 
conv <— 1 psi 
j^O 

while j < 20 
I for i € 0.. 2 


-P. R 


a <Wso 4 ---[l + r iH^ + v H £i"| 

l P 

-P.R 

o«si<-;-[l +r iHm- v-He'1 

t n -* 


-P i R 

r xxso 4 -*0.3 + r H £» 


-PR 

a xxsi 4 -‘0.5 - T-H £» 


r <Wfo ' 


a Wtfi <“ 


-PR 

l 


-P. R 

l 


i - r 


l - r 



3 \ 05 71 

1 + V 

—K 


, 2 

- 

U-v J JJ 


1 - V 


CT xxfo 


-P-R 

i 


1 - V 

05 


05 


K 


a xxfi 4 “ 


-P. R 

l 


o.5 - r 


o.5 + r 


U-v 


1 - V 


K 


,0.5 


■K 


a gy +- 


^ ® (#60 ^ 


a #si 
a xxso 
a xxsi 
a #fo 
° #fi 
a xxfo 
V a xxfi; 


a 2 <-o sy 


a 3 «- O' sy, 


17 4 <- o sy 
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Perl ~ 

Perit : 


a SYM ' 


f 2 A 

\° 1 -oia 2 + o 2 J 


( 2 2 ^ 

o SYF <- \o 3 - 03 - 04+04 ) 

f symV 

SY <- max 

[\^SYF ) 

stress. <— a sy 

test. <— (f y - a SY! 

test 0 


shel <— 


dav <- 


conv 
testj 

conv 


break if jstressj - fy| < limit 
?0 ^ shel-dav > 0 

?2 <— P| if shel-dav < 0 


P 1 *- p o + 


i ^2 P 0 - 


j<- j+ 1 


OUtg 4- 

1 


1 1 


conv 
out* <— j - 1 


out 


^ 819.180 
V 14 J 

Pcrl o -lpsi 
Pent = 819.2psi 
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PART 2 LOBAR BUCKLING 


Windenberg Approx of Von Mises (1933) 

Assumes n lobes = Pi*D/L 

Collapse 

pressure: 



p cLB = 2947.7psi 

PART 3 GENERAL INSTABILITY 

Corrected Bryant Formula (1954) for better model test correlation 
Pressure loading is: 

P:= p g D t SF g j 
P= 1786.13 lpsi 

Compute effective frame spacing: 



y = 0.473 

Compute clear length: 
Lc := Lf - t w 
n| := 0.5-^ 1 - y 
n] =0.363 
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n 2 := 0.51 + y 
ri2 = 0.607 

2 2 
cosh^nj-O) -cos(n2-0) 

cosh (nj-0)sin h(ni-e) cos(n2-0)-sin(n2-0) 

n l n 2 

Fj = 0.894 must be less than 1.00 

F) is almost a linear decreasing function for pressures from 1 to 2000 psi with an average value 

of 0.90. This will be used in the following analysis as the pressure is the unknown and therefore 

the above equations cannot be directly used. 

F] := 0.90 

Leff := L C F] + t w Effective shell plate length: 

Lgff = 1.234in 

Theoretical critical lobe number values are: i := o.. 2 
A e ff := L e fftp Effective plate area 

2 

Aeff “ 0.106in ^ 

Circumferential Lobes 



Longitudinal Lobes: 
R 

m:= n — 

4 

m= 3.062 


4 
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I e ff = 7.711x 10 _3 in 4 

Determine Rf as radius to centroid of combined plate and stiffener: 
The critical Elastic General Instability pressure is: 
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Min Pressure: 


PcGl 


4 \ 

2.039x 10 
8.546x 10 3 
^6.379x 10 3 , 


psi 


f 2^ 


n := 


3 

< 4 > 


PcGI := min(P cGI ) 


P C GI = 6378.9psi 


Summary 

Axisymmetric Buckling 
Peril = 819.2psi 
Lobar Buckling 
p cLB = 2947.7psi 
General Instability 
P c Gi = 6378.9psi 
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Analytic Solution 


Defintions 

ksi := 6.89475710 6 Pa 
E := 3000(ksi 
R := 8.4179n 

D := 2 R 
Oy := 65.5ksi 
L:= 1.36fin 
Lj, := 8.636n 

p:=.3 

t := 0.0858n 
Ring Stiffeners 

t w := 0.044n 
H := 0.5319n 
b:=t w 
bf:= 0.399n 
b 2 :=bf-t w 
d := 0.078n 
4)nash := L - bf 


Modulus of Elasticity 

Radius of cylinder to centerline of shell 

Diameter of cylinder 

Yield strength 

Length of supported cylinder 
Distance between bulkheads 

Poison's ratio for Fe/Steel 
Shell thickness 

thickness of web of ring stiffener 
height of ring stiffener 

faying width of stiffener (from P&S for I beam stiffener) 
breadth of ring stiffener 

breadth of ring stiffener minus the web thickness 
ring stiffener flange thickness 
distance from flange edges 
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Cl := 5 


V H+b r d 


first centroidal height of ring stiffener 


C] = 0.143in 


c 2 :=H-Cj 


second centroidal height of ring stiffener 


C2 = 0.389in 


h := cj - d distance from centroid of ring stiffener to nearest edge of flange 


h = 0.065in 


Rj. := R + ,5-t + c 2 radius to centroid of ring stiffener 


Rj. = 8.85in 


A r := (t w H + b 2 d) cross-sectional area of ring stiffener 


A r = 0.051 in 2 


V := j-C bf * c l 3 “ b 2* h3 + 


moment of inertia of ring stiffener about its centroidal axis 


l T = 1.219x 10 3 in 4 


^eff - ^r’ 


effective area of stiffener eqn [24a] from P&S 


A e ff = 0.049in 


ratio of effective frame area to shell area eqn [62] P&S 


a =0.415 




ratio of faying width to frame spacing eqn [62] P&S 


p = 0.032 
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Combined Plate and Ring Stiffener 

A p := L t area of effective plate 

A sp ;= A p + A r area P^ ate and ring stiffener 

H c := H +1 height of combined plate and ring stiffener 

H c = 0.051ft 

B ] := L - t w plate length minus the web thickness 

t w -H c 2 + B, t 2 + b 2 d (2 H c - d) 

c ic := 5 “ +Bi . 1 + d . b - neutral axis of combined plate and ring stiffener 

from outer fiber of plate (R 0 ) 

Cj c = 0.174in 

R c := R - ,5-t + c lc radius to centroidal axis of combined ring stiffener and shell 

Rj. = 8.549in 

T e := 3 '[ Lc 1 c 3 " B l ( c lc - l ) 3 + bf ( H c - c lc) 3 - ( bf - ^ K - c lc ~ d f 
moment of inertia of combined plate and shell 
l e = 0.0 liif* 

a. Buckling of unreinforced shells 


guess at number of waves around the circumference 
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Given 
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b. Axisymmetric Buckling 


Pcab -= 


y 4- 0 

limits- lOpsi 
test 4- Opsi 
conv 4- lpsi 
j 4 0 

while j < 20 

n 2 ^ 


n — 

V3U-, 


2 L 


yfR -t 


Fi <- 


cosh ^ |-0j 

| 2 -cos( 

n 2 e) 2 

cosh^ri j 0j 

|-sinh| 

(n | 0 ) cos 

M 

I’H 

hr 0 ) 


cosh 


1| 1 2 
(n ] g) sin^n2 °) sinh(n j ©) coso) 


^2 




cosh^B j- 0 j sinh(n | 0 ) cos2 ^)* sin (H2'®) 
ill n 2 

cos^n2 g )‘ sin (n2* g ) coshidj-sinh^ j-oj 

3 n 2 0i 


Fj- P 

>( I- 


3 

\>V 


^2 cosh j-Oj-sinh^r) j-0 j cos(n 2 'G) sin(n 2 oj 

ni n 2 

cosh^n j'Oj-sin^r) 2 o) sinh^n i'6)cos(r| 2 -0) 

12 n I 


cosh(r| j-0j sinh(n |-0j cos (l2' 0 ) •sin(n 2 e) 

1 | n 2 


Hh 

€» + (J + (l - p)-F, 


denoml 4 — A 


denom2 < 


F 2 2 + F 2 F 4 (l-2 M ) 


( 

091 1 

1 t 2 f 

2 ) 

U 

.-m 2 J 

+ F 4 *\ 

* - M + M /• 


0.91 


1 - M 



Pc2 * 


•a|f 2 ~ m .f 4 


0.91 ' 

J 1 -n 2 . 


3 

N 4 


+ denoml - denom2 
break if |p c2 - test | < limit 

P C 2 


2-H 
test 4— p c2 

j4- j+ 1 

Pc2 


|[jnr7)](i) 


out„ 


conv 


out 4 - j 
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d. General Instability of shells and rings 



Given 

n > 1 
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n 2 := Minimiz^p c Qj,n) 

n 2 = 3.957 n must be an integer 

n 2int :=round ( n 2>°) 


n 2int = 4 


1 2prime* 


n 2int “ 1 
n 2int 
n 2int + 1 
n 2int + 2 


p cgi2 := 


PcGl^ n 2prim^ 
PcG]^ n 2primej^ 
PcGI^^prime,^ 
PcGl^ n 2prime^ 

Pcgibryant -= m ' r (Pcgi2) 


Pcgibryant = 6498x 10 P si 


Summary 

Lobar Buckling 
PcLB = 2947.7p si 
Axisymmetric Buckling 
Paxiy = 8? 5 9p si 
General Instability 
Pcgibryant = 6498.2psi 
n 2int = 4 
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Appendix D: Codes for Test Cylinder 2.c 
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AMERICAN BUREAU OF SHIPPING 


Rules for Building and Classing Underwater Vehicles, Systems, and Hyperbaric Facilities 

Definitions 


ksi := lOOC^si 
E := 3000(ksi 

v := 0.3 

Oy := 15700Q)si 
Shell Parameters 
L s := 3.256n 
L c := 115.533n 


R := 18.882n 


t := 0.337ln 


Modulus of Elasticity 
Poisson's Ratio 
Yield strength 

Distance between stiffeners 
Distance between bulkheads 

Mean radius of shell 
Thickness of shell 


R o : =R + 


t 

2 


Outer radius of shell 


R 0 = 19.051 in 
D o ; = 2 *o 


Outer Diameter of shell 


Ring Stiffeners 

t w := 0.122n 


depth := 2.315n 


b 


t 


vv 


bf := 1.552n 


b 2 := bf - t 


vv 


thickness of web of ring stiffener 
height of ring stiffener 

faying width of stiffener (from P&S for I beam stiffener) 
breadth of ring stiffener 

breadth of ring stiffener minus the web thickness 
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d := 0.305n 


ring stiffener flange thickness 


:= l s - t w Unsupported spacing between stiffeners 

Ljj = 3.129in 
L:= ma^L^Lj,) 

L = 3.256in 

( 2 

t w -depth + b 2 'd 

ci := .5- - first centroidal height of ring stiffener 

^ t w - depth + b 2 -d J 

Cj = 0.558in 

c 2 := depth - cj second centroidal height of ring stiffener 

c 2 = 1.757in 

h := cj - d distance from centroid of ring stiffener to nearest edge of flange 

1*3 := r - ,5-t - c 2 radius to centroid of ring stiffener 
Rg = 16.957in 

A s := (t w depth + b 2 -d) cross-sectional area of ring stiffener 
A s = 0.729in 2 

Ij. := ^j-^bf-C] 3 - b 2 -h 3 + t w c 2 3 j moment of inertia of ring stiffener about its centroidal axis 
I r = 0.312in 4 

Rf := R - ,5-t - depth Radius to tip of the stiffener 

Rf = 16.398in 
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Combined Plate and Ring Stiffener 


effective plate length 

L e = 2.442in 

A p := L t area of effective plate 

A sp := Ap + A s area of plate and ring stiffener 

H c := depth + t height of combined plate and ring stiffener 

H c = 0.221 ft 

B| := L - t w plate length minus the web thickness 

t vv H c + B r l + ~ 

c ic := 5 --—--—- neutral axis of combined plate and ring stiffener 

t w H c + Bjt + db 2 



C| c = 0.937in 


from outer fiber of plate (Ro) 


Rc ;= R + -5 1 - c Ic radius to centroidal axis of combined ring stiffener and shell 


R c = 18.114in 

! e := J [ L c lc 3 - B l ( c lc " l ) 3 + bf ( H c " c lc) 3 " ( bf - * w HHc - c lc - d) 3 ] 
moment of inertia of combined plate and shell 
I e = 1.946in 4 

General Equations 



M = 1.291 
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0 := 



■M 


0 = 1.659 


Q:= 


0 _ 

2 


Q = 0.83 

cosh (2-Q) - cos(2 Q) 
sinh(2-Q) + sin(2Q) 

N = 0.797 


G . = 2 (sinh(Q) cos(Q) + cosh(Q) sin(Q)) 
sinh(2-Q) + sin(2*Q) 

G = 0.926 

^ sinh(2Q) - sin(2*Q) 
sinh(2-Q) + sin(2Q) 


H = 0.435 


Inter-Stiffener Strength (6.19.1) 

1) Inter-stiffener strength equations 

yield pressure at midbay and midplane of cylinder 

This equates to axisymmetric buckling 


A := A s‘ 


R 


Effective area of plate and stiffener (Internal stiffeners) 


A = 0.81 tin 2 
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V_2 


F := ■ 


A + t w -t + 


2Nt-L 

6 


F = 0.335 


P ./A 

y 1 - F 


P y = 421 l.lpsi 


von Mises buckling pressure for a cylinder 
This corresponds to Lobar buckling 


P 


m 


5 



P m = 13410.8psi 



if 1 < 


m 


< 3 




P c = 3509.3psi 


maximum allowable working pressure for inter-stiffener strength 


p aits ■ p c - 8 


P aiK = 2807.4psi 


2) longitudinal stress 


limiting pressure corresponding to the longitudinal stress at stiffeners reaching yield 
No direct correlation to major failure modes 


A + hv' 1 + 


2-N-t-L 


ZOy-' 12 V 

-—• 1 + - yH 

R , 2 
U-v / 


Pi = 3566. lpsi 


maximum allowable working pressure for longitudinal stress 


p als : = p l- 67 


p a]s = 2389.29psi 


b Overall Buckling Strength (6.19.5) 


A 2 (n) := n - 1 
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Et E1 e A 2( n ) 

Pnl( n ) ;= V A l (n) + -;- 

R r 3 l 

Given 
n > 2 

n 2 •= Minimiz^p n | ,n) 

n 2 = 2 n must be an integer 


n 2int :=round ( n 2’°) 
n 2int = ^ 



Pn :=min (Pn2) 

p n = 8642. lpsi 
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Summary 

Axisymmetric Buckling 
P y = 421 l.lpsi 
Lobar Buckling 
P m = 13410.8psi 
General Instability 
p n = 8642. lpsi 

n 2int = 2 




Germanischer Lloyd 

1988 Edition 

Chapter 2 - Submersibles 

MPa := 10 6 Pa 


ksi := 10 


psi 


General Definitions 


E:= 3000Cksi 
Lj := 3.256n 


Modulus of Elasticity 
stiffener spacing 


Lj := 115.532n 


s := 0.337in 


R := 18.882n 


v := 0.3 

k := 157 ksi 


length of cylinder between bulkheads or lines of support 

thickness of shell 

radius to centerline of shell 

Poisson's Ratio 

minimum yield stress of material 


Stiffener Dimensions 


0.305n 
d f := 0.7125n 

b2 

d w := 2.01in 
e w := 0.121n 


flange thickness 

width of flange from web to edge of flange 

Height of web 
web thickness 

width of stiffener ring in contact with shell 


192 




Aj := e f (e w + 2d f ) + e w -d w cross-sectional area of stiffener ring 

Aj=0.729in 2 

H stiff := d w + e f 

H stjff = 2.31 Sin 


1 e w H stiff +b 2 e f 


1 2 e w' H stiff + b 2 e f 


C| = 0.558in 


c 2 H stiff “ C 1 


C 2 = 1.757in 


h 2 C 1 _ e f 


h 2 = 0.253in 


*1 : = 


-Y ^e w + 2 d f ) cj 3 - (2-d f ) h 2 3 + e w -c 2 3 J Centroidal Moment of inertia of ring stiffener 


= 0.312in 


e := c 2 + .5 s distance from stiffener centroid to center of shell 


e = 1.925in 

Effective Stiffener and Shell 

L etest : = b + 


L etest = 93.838mm 


V= 


L etest lf L etest - L 1 


Effective length of shell 


Lj otherwise 


L e = 3.256in 
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I e = 1.946in 


Moment of Inertia of combined plate and shell 


^stiffc • H stiff + s 


Total height of stiffener and plate 


H stiffc = 2652in 


d,:=s 


bs ,:= L e-e w 


2 2 

1 e w'^stiffc + b 2 e f + bs r s ( 2 '^stiffc ~ s ) 

2 e w' H stiffc + b 2 e f + bs ] s 


C 1 C ; =- 


Cj c = 1.715ln 


c 2c' ^stiffc c lc 


c 2c = 0.937in 


Rq := R + 5 s - c 2c radius of stiffener ring centroid including effect of L e 


R 0 = 18.114in 


A := Aj - 


Modified area of stiffener ring 


A = 0.792in 2 


e, :=c 2c + ,5 s 


ej = 1.105in 
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2.4 Asymmetric Buckling (Lobar but not named this) 



n := 2 initial guess 

s 2 (n 2 - 1 + A.j 2 ) 

12 R 2 (l - v 2 ) 

’ ' 2,2 
n - 1 + .5 X X 

Es-P n] (n) 

Pnl( n ) := --- 

Given 



n > 2 

n g :=Minimiz^p n ],n) 

n g = 10.171 n must be an integer 

n gint := round (n g ,o) 


n gint 


10 


n gfirst • 


n gint 

n gint 


100 if n gint= 1 
- 1 otherwise 


n gprime 


n gfirst 
n gint 
n gint + * 
n gint + 2 




<P nl( n gpi-ime 0 ) 

p nl^ n gprimej j 

p n2 := / \ 

Pnl^gprimejJ 

< Pn l( n gpnme 3 j 


p cr :=min ( p n2) 


p cr = 12139.3psi 


3.3.3 General Instability 


L:=L3 


Spacing between effective stiffeners 


Pn2<"):= 


n — 1 + 0.5 -^2 J*l n + 


2 2 .2 


P 0 (") - 


Esp n 2( n ) 


Pl(n) := 


(n 2 - 1)-E-1 


R o L 1 


p nla( n ) := p 0 ( n ) + p l( n ) 


n 2 :=Minimiz^p n]a ,n) 


n 2 = 2 n must be an integer 
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1.768x 10 5 
9.702x 10 3 
2.418x 10 4 
4.525 x 10 4 


psi 


2.6 Symmetric Buckling (axisymmetric) 

More Definitions 

For Elastic-Plastic Region (Eqns are in the program) 

E s := E Secant Modulus = Young’s Modulus for elastic region 

Ej := E Tangent Modulus = Young's Modulus for elastic region 

v n := v Poisson’s Ratio, elastic-plastic; = Poisson's ratio in elastic 

shape factor 

a = 0.02— 
mm 
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F 4 := 


cosh (Cg) sin sinh (Cg) ■ cos (c^ 

3 C 7 _ C 6 

j 2 cosh(Cg)-sinh(Cg) cos(c^-sin(C9) 


c io : -‘ 


1 - - 
2/ 


s-Lj 


A b ( b ^ 
+ — + | 1 - 

L iy 


s-Lj Lj 


F 1 


c n 


0.91 


1 -v 


a o := ' 


-p-R 


a x :_a o' - + C 10 C 11 F 4 


V^o’O " C 10 F 2 + v C 10' C H F 4) 


c xl : - CT o’ 2 _C 10' C 11' F 4 


a x = -1.168x 10 Pa 




1 _ c io f 2 _ vC io c i r F 4; 


a ( j ) = ”1.259x 10 Pa 


K 0 := ' 


2 2 

G i : ’J G x + a <|> " a x ff 4 


a- = 176.383ksi 


a:=-K 0 a x 


a = 182.6ksi 
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^ = 1 . 067 x lo'°Pa 
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trial := 


a r L i 

7t 

trial = 0.413 

iterate on m until trial <= trial2 


m := 1 


ria!2:= 


m 


tria!2:= f —*(m+ 1) 
2 


trial2= 1 


Pm* Paa‘^0' 


/ ct r L i A 

V Tt-m ) 


1 


r 7i-m > 

K a r l \j 


p = 4566.9psi 


p s 456^)si 


iterate on p until p m equals p 


Summary 

Axisymmetric Buckling 
p m = 4566.9psi 
Lobar Buckling 
p cr = 12139.3psi 
General Instability 
Pg = 9702.4psi 
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Society of Naval Architects and Marine Engineers 


Submersible Vehicle Systems Design 
ksi := lOOQjsi 

Input Section 


Lf := 3.25fin 

Length between frames 

:= 115.532n 

Length between bulkheads 

R := 18.882n 

Radius of cylinder to centerline of shell 

t := 0.331n 

shell thickness 

Rq i= R + .5-t 

radius to outside of shell 

Rq - 19.05 lin 

E:= 30000ksi 

modulus of elasticity 

Gy := 157 ksi 

minimum yield stress of material 

\i ~ 0.3 

Poisson's ratio 

Ring Stiffeners 

t w := 0.1 21n 

thickness of web of ring stiffener 

depth :=2.315n 

total height of ring stiffener 

bf := 1.552n 

breadth of ring stiffener 

d := 0.305n 

ring stiffener flange thickness 


b 2 := bf - t w 


breadth of ring stiffener minus the web thickness 



first centroidal height of ring stiffener 


( 2 2 
t w -depth + b 2 d 

c,:=.5- - 

^ t w -depth + b 2 d 

Cj = 0.558in 

c 2 := depth - cj second centroidal height of ring stiffener 

C 2 = 1.757in 

h := c | - d distance from centroid of ring stiffener to nearest edge of flange 

h = 0.253in 

R,. := R - .5't - c 2 radius to centroid of ring stiffener 

R,.= 16.957in 

A r :=(t w depth + b 2 -d) cross-sectional area of ring stiffener 

A r = 0.729in 2 

3 3 3 ^ 

bf c, - b 2 h + t w c 2 J moment of inertia of ring stiffener about its centroidal axis 
l r =0.312in 4 

Combined Plate and Ring Stiffener 

A p := Lpt area of effective plate 

A S p := A p + A r area °f P^ ate anc * ring stiffener 

H c := depth + t height of combined plate and ring stiffener 

Bj := Lf - t vv plate length minus the web thickness 
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neutral axis of combined plate and ring stiffener 
from outer fiber of plate (Ro) 

Cj c = 0.937in 

r c := r + ,5-t - c lc radius to centroidal axis of combined ring stiffener and shell 

R,, = 18.114in 

'e : =| [ L F<=lc 3 - B l ( c Ic - •)’ + bf ( H c - c lc) 3 - ( bf - >wM h c “ c lc ' “f] 
moment of inertia of combined plate and shell 
I e = 1,946in 4 


a. Axisymmetric yielding 



0 = 1.659 

cosh(0) - cos ( 9 ) 
sinh(0) + sin(9) 

N = 0.797 


V H c + B l-t + b 2 d ( 2 ~ H c ~ d ) 
V H c + B r t + db 2 
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py = 3863.5 Ipsi 


b. lobar buckling 



Pj, = 1.341 x 10%si 


c. General Instability 

7t-R 
m:=- 

Lb 

n := 1 
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PcrlO) 



Given 

n > 1 

n 2 := Minimiz^p cr |, n) 

n 2 = 1 -654 n must be an integer 

n 2 int :=round ( n 2 ’°) 
n 2 int = 2 


n 2 first : “ 


n 2 in t ] 0 ° if n 2im = 1 


n 2 int “ * otherwise 


n 2prime’ 


n 2 first 
n 2 int 
n 2 int + 1 
v n 2 int + 2 ) 





J 


p cr = 8642. lpsi 
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Summary 

Axisymmetric Buckling 
Py = 3863.5psi 
Lobar Buckling 


p b = 13410.8psi 


General Instability 
p cr = 8642. lpsi 
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MIT 13A Professional Summer Submarine Design Trends 

PROGRAM TO COMPUTE SUITABILITY OF SUBMARINE DESIGN PARAMETERS. 

Ref: "Hull Material Trade Off Study", D Fox, Jan 94 
Define input parameters: 

ksi = 100Q)si 

Global Variable Inputs: 


e = 0.0 in 

Eccentricity 

Material: 


Gy = 15700Q)si 

Yield Strength 

E = 3000Cksi 

Modulus of Elasiticity 

to 

o 

III 

> 

Poisson’s Ratio 

Geometry: 


R = 18.882n 

shell radius 

Dh2-R 

shell diameter 

Lf = 3.256n 

frame spacing 

4 = 115.532n 

bulkhead spacing 

tp = 0.33^n 

shell thickness 

tf = 0.305n 

flange tickness 

Wf = 1.5 525n 

flange width 

t w s 0 . 12 ^n 

web thickness 

h w = 2 . 01 in 

web height 


209 



Compute areas: 

1 p 

Ri := R - — Frame Radius 

2 

R] = 18.713in 

Frame flange, web area: 

Af := tf-Wf 

A\V := Uv’^w 

A := Af + a w Frame Area = Flange + Web 

A = 0.729in 2 


PARTI SHELL YIELDING 

bv'tp 

B := - Area ratio 

A + t w -tp 

B = 0.055 



cosh(o) - cos(6) 
sinh(9) + sin(o) 


Slenderness parameter: 


Deflection coefficient: 


p:= 


2-N 

A + t w -tp 




Frame flexibility parameter: 


P = 1.366 
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PART 3 GENERAL INSTABILITY 

Corrected Bryant Formula (1954) for better model test correlation 
Pressure loading is: 

P := p-g-DfSFgj 

P= 1.674x 10 3 psi 

Compute effective frame spacing: 



y = 0.145 
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Compute clear length: 
he := Lf - t w 
nj 0.5*>/1 - y 
ni =0.462 
Web thickness: 
n 2 := 0.5^1+ y 
n 2 = 0.535 
t w = 0.127in 

2 2 
cosh^nj-0^ - cos (n 2 0) 

cosh(n|-0)sinh(n]-0) cos (n 2 -0) sin(n2-0) 

n l n 2 

F| = 0.959 must be less than 1.00 

FI is almost a linear decreasing function for pressures from 1 to 10000 psi with an average value 
of 0.96. This will be used in the following analysis as the pressure is the unknown and therefore 
the above equations cannot be directly used. 

F! := 0.96 

Effective shell plate length: 

Left := Lc'Fl + »w 

L e ff=3.131in 

Theoretical critical lobe number values are: i : = o.. 2 
A e ff := Lgff t p Effective plate area: 

A e ff= 1.055in^ 
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Circumferential Lobes 
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Total: 


Jeff := Ipcor + Kvcor + Ifcor Moment of Inertia corrected for neutral axis. 


leff = 1.92 in 4 

The critical Elastic General Instability pressure is: 



Min Pressure: 



( 3> 




8.536x 10 


f2\ 

p cGl = 

2.107x 10 4 

psi n := 

3 


^3.942x \0 4 j 


u, 


P C GI := min(P cGI ) 

P C GI = 8535.8psi 
n gi := 2 

Summary 

Axisymmetric Buckling 
Pcrit = 37I2psi 
Lobar Buckling 
p cLB= 13410.8psi 
General Instability 
p cGI = 8535.8psi 
n g i= 2 
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Analytic Solution 


Definitions 

ksi := 6.89475710 6 Pa 
E := 3000(ksi 

R := 18.882n 

D := 2R 
Gy := 157ksi 
L := 3.256n 

:= 115.532n 

M:= -3 

t := 0.332n 

Ring Stiffeners 

t w := 0.121n 
H := 2.315n 
b: = t w 
bf := 1.552n 

b 2 ;= bf " l w 
d := 0.305n 

Lbnash ;= L _ bf 


Modulus of Elasticity 

Radius of cylinder to centerline of shell 

Diameter of cylinder 

Yield strength 

Length of supported cylinder 
Distance between bulkheads 

Poison's ratio for Fe/Steel 
Shell thickness 

thickness of web of ring stiffener 
height of ring stiffener 

faying width of stiffener (from P&S for I beam stiffener) 
breadth of ring stiffener 

breadth of ring stiffener minus the web thickness 

ring stiffener flange thickness 
distance from flange edges 
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c,:= .5- 


V H 2 + b 2 -d 2 
t w 'H + b2*d 


first centroidal height of ring stiffener 


Cj = 0.558in 


c 2 := H - cj second centroidal height of ring stiffener 


C 2 = 1.757in 


h := Ci - d 


distance from centroid of ring stiffener to nearest edge of flange 


h = 0.253in 


R,. := R - 5 t - c 2 radius to centroid of ring stiffener 


R,. = 16.957in 


A r := (t w H + b 2 d) cross-sectional area of ring stiffener 


A r = 0.729in 2 


f 3 ) ( bf C l 3 ~ b2 ' h3 + l w c 2 3 ) moment of inertia of ring stiffener about its centroidal axis 


I r = 0.312in 4 


A eff := A r 


effective area of stiffener eqn [24a] from P&S 


Aeff = 


ratio of effective frame area to shell area eqn [62] P&S 


a = 0.739 


P :=- ratio of faying width to frame spacing eqn [62] P&S 


P = 0.039 
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Combined Plate and Ring Stiffener 


A p := L t area of effective plate 

A := A + A area of plate and ring stiffener 

sp p r 

H c := H + t height of combined plate and ring stiffener 

H c = 0.221ft 

Bj := L - t w plate length minus the web thickness 


c Jc := .51 


V H c 2 + B r t2 + b 2 d ( 2H c- d ) 


l w' H c + B l l + db 2 


neutral axis of combined plate and ring stiffener 
from outer fiber of plate (Rq) 


c lc = 


R C :=R+ .5t-c lc 


radius to centroidal axis of combined ring stiffener and shell 


R c = 18.114in 

h : = j[ L c lc 3 - B l ( c lc - O' + bf ( H c - c lc) 3 " ( bf - *wK H c - c lc - d ) 3 ] 

moment of inertia of combined plate and shell 


I e = 1.946in 4 


a. Buckling of unreinforced shells 

von Mises bucking pressure: 

n := l guess at number of waves around the circumference 
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Purl< n ) := 



Given 


n > 1 

n i := Minimiz^p ur j ,n) 

n] = 10.14 n must be an integer 

n lint :=round ( n l>°) 



4 

p ur = 1.214x 10 psi 
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b. Axisymmetric Buckling 
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f 4080.3^) 

Pcab = - 

\ 2 ) 

Paxiy := Pcab Q 'P si 

Paxiy = 4080 - 3 P si 

c. Asymmetric Collapse (Lobar buckling) 

Windenberg Approx of Von Mises (1933) 
Assumes n lobes = Pi*D/L 
Collapse pressure: 




Given 

n > 1 
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n 2 := Minimiz^p c Qj,n) 

n 2 = l .654 n must be an integer 


n 2int : = round ( n 2’°) 


™2int : 


1 2 prime‘ 


^ n 2int - 1 ^ 
n 2int 
n 2int + 1 
V n 2int +2 ; 


Pcsi2 • 


PcGI^ n 2prime^ 
PcGl^prime,) 
PcGI^ n 2prim^j 
PcGI^ n 2prime^ 

Pcgibryant := mi,1 (Pcgi2) 


Pcgibryant = 8 ' 642x 10 P si 


Summary 

Lobar Buckling 
PcLB = 13410.8psi 
Axisymmetric Buckling 

Paxiy= 4080 - 3 P si 

General Instability 

Pcgibryant = 8642 1 P si 
n 2int = 2 
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